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Highlight  nummary 

This  report  gives  a  complete  nummary  of  the  tent  data  on  Hhaler. 
generated  by  ISHC  and  Purdue  University.   Details  are  Riven  for  a  very 
simple  storage  and  retrieval  system  which  is  adequate  for  the  small 
amount  of  data  (l63  sets)  available  at  the  present  time. 

Results  are  presented  from  the  different  statistical  analyses  that 
were  performed.   These  results  include  histograms ,  bivariate  correlation 
coefficients  and  regression  equations.   Reasonably  good  bivariate  corre- 
lations exist  between  the  different  indices  describing  the  slaking  resis- 
tance of  shales.   These  correlations  are  however  improved  by  using  quadratic 
equations.   Various  regression  equations  are  proposed  for  determining  CBR 
from  various  parameters,  usually  a  combination  of  five. 

It  became  clear  during  the  investigation  that  it  is  important  to 
have  as  many  complete  data  sets  as  possible  for  future  analyses.   The 
standardization  of  testing  methods  is  also  of  utmost  importance  in  order 
to  increase  the  potential  of  the  data  bank. 


Laboratory  Test  Results  on  Indiana  Shales: 
Storage,  Retrieval  and  Statistical  Analyses 

1 .   Introduction 

Large  amounts  of  soil  test  data  are  being  generated  in  every  state 
during  the  testing  of  material  for  new  highway  facilities.   Centerline 
samples  and  samples  of  proposed  material  sources  are  subjected  to  series 
of  'standard'  material  tests.   The  test  data  usually  find  their  way  to 
the  archives  or  the  waste  paper  basket  after  completion  of  the  design 
and  construction  of  a  facility.   The  usefulness  of  these  generated  data 
was  not  realized  until  recently  and  a  number  of  states  are  now  in  the 
process  of  introducing  storage  and  retrieval  systems  (Crawford  et  al . , 
1972;  Spradling,  1976).   Such  a  system  is  easily  handled  by  a  computer, 
especially  when  a  large  number  of  data  sets  are  available.   For  a  storage 
and  retrieval  system  to  be  as  complete  as  possible  it  is  necessary  to  in- 
clude data  generated  by  everybody  involved  in  material  testing,  e.g., 
state  highway  materials  laboratories ,  consulting  engineers  and  research 
organizations . 

It  is  very  important  that  a  storage  and  retrieval  system  be  central- 
ized so  that  potential  users  can  obtain  the  data  in  which  they  might  be 
interested.   The  availability  of  such  a  data  bank  should  be  brought  to 
the  attention  of  potential  users. 

The  type  of  storage  and  retrieval  system  to  be  used  depends  on  the 
amount  of  data  available.   Any  system  should  be  flexible  enough  so  that 
changes  can  be  made  at  a  later  stage  with  the  least  amount  of  effort. 
Any  storage  and  retrieval  system  should  therefore  be  seen  as  a  dynamic 


system  because  it  will  be  changed  from  time  to  time  to  fit  the  needs  of 
the  users  and  to  make  its  administration  easier  and  more  economical. 
These  qualities  can  be  best  incorporated  in  a  computerized  system. 

The  statistical  analysis  of  stored  data  can  lead  to  a  number  of 
very  useful  results.   Experience  usually  indicates  some  'typical  proper- 
ties and  behavior'  of  certain  types  of  material.  Statistical  analyses 
can  prove  or  disprove  such  notions.   Statistical  analyses  will  also 
indicate  the  typical  range  of  values  to  be  expected  for  different  para- 
meters of  the  different  members  in  the  system  as  well  as  the  distribution 
of  these  parametero-  This  can  play  a  very  important  role  in  assessing 
the  suitability  of  a  certain  sample  before  elaborate  laboratory  testing 
is  undertaken.  Statistical  analyses  should  be  updated  (regularly)  as 
more  data  become  available.  They  are  therefore  dynamic  and  the  results 
can  be  stated  with  graater  confidence  with  increase  in  the  size  of  the 
available  data  bank.  The  results  are  just  as  good  and  complete  as  the 
data  available  at  the  time  of  the  analyses.  An  additional  benefit  is 
the  development  of  new  ideas  for  more  productive  research  based  on  em- 
pirical relationships  and  correlation  of  soil  properties  (Spradling,  1976) 

Correlations  between  different  parameters  can  be  developed  by  using 
stored  data.  Regression  models  can  be  established  to  predict  the  para- 
meters which  are  difficult  to  measure  by  using  a  set  of  simple  parameters. 
These  models  should  be  evaluated  as  new  data  become  available  to  test 
their  reliability. 

This  particular  study  was  conducted  to  set  up  a  storage  and  retrieval 
system  for  the  test  data  generated  on  Indiana  shales.  This  system  was  set 
up  separately  from  a  comprehensive  storage  and  retrieval  system  being 


introduced  presently  by  the  Indiana  State  Highway  Commission.   The 
purpose  of  this  system  was  to  be  of  help  to  the  Indiana  State  Highway 
Commission  in  providing  them  with  information  on  shales  while  they  are 
setting  up  their  own  storage  and  retrieval  system.   This  system  will 
also  be  an  aid  in  the  development  of  research  projects  on  shales  at 
Purdue  University.   It  is  hoped  that  useful  practical  results  will  also 
be  obtained  from  the  statistical  analyses  done  during  the  course  of 
this  study. 

Certain  descriptors  had  to  be  used  in  order  to  locate  the  different 
data  sets.  In  the  system  described  in  this  report,  laboratory  number, 
geographic  (counties),  geological  (system,  series,  stage)  and  physio- 
graphic (bedrock  units)  descriptors  are  used.  Retrieval  should  be  done 
in  the  most  useful  way,  e.g.,  based  on  alphabetical  order  for  certain 
parameters  and  numerical  order  for  others.  This  type  of  ordering  will 
help  to  increase  the  usefulness  of  the  system. 

The  sources  of  data  for  the  system  described  in  this  report  are  the 
testing  done  by  the  Indiana  State  Highway  Commission  on  shales  as  well 
as  the  data  given  in  the  three  theses  dealing  with  shales  submitted  to  Purdue 
University  (Deo,  1972;  Chapman,  1975;  Bailey,  1976). 

This  report  describes  the  information  stored  and  the  retrieval  of 
the  data.  The  results  of  various  statistical  analyses  are  also  presented. 
Finally,  recommendations  are  made  for  future  research. 

2.   Geology  of  Indiana  Shales 

The  information  of  this  section  was  taken  from  the  thesis  by  Deo 
(1972).  Another  useful  reference  in  this  respect  is  that  by  Harrison  and 
Murray  (196U).  A  geological  map  of  Indiana  is  given  in  Figure  1. 
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Figure  1.   Bedrock  geology  of  Indiana  (after  Wayne,  1956). 


In  many  parts  of  Indiana,  shales  are  either  exposed  at  the  earth's 
surface  or  underlie  it  at  shallow  depths  that  are  within  the  range  of 
engineering  considerations.   Only  shales  of  the  Palezoic  Era  are  present 
in  Indiana  and  hence  the  montmorillonitic  clays  related  to  more  recent 
rocks,  volcanic  activity,  and  weathering  in  arid  regions  are  not 
represented . 

2.1  Shales  of  Ordovician  Age 

The  oldest  geologic  system  of  rocks  in  Indiana  that  contain 
shale  of  engineering  significance  is  the  Ordovician.  These  rocks  are 
exposed  in  the  southeastern  part  of  the  state.   One  such  rock  unit  of 
importance,  the  "Dillsboro  Formation",  lies  within  the  Ordovician. 

a)  Shale  of  Dillsboro  Formation.   This  formation  consists  of 
alternating  beds  of  shale  and  limestone;  at  some  locations  more  than 
five  hundred  beds  of  alternating  shale  and  limestone  can  be  observed." 
The  thickness  of  shale  beds  varies  between  one  inch  and  two  feet. 

The  Dillsboro  Formation  lies  above  the  Kope  Formation  and  directly 
below  the  Saluda  Formation.   Limestone  in  the  Dillsboro  is  argillaceous 
and  shales  are  calcareous.   Common  clay  minerals  present  are  illite, 
kaolinite  and  chlorite. 

Generally  the  shales  are  highly  fissile,  and  with  repeated  wetting 
and  drying,  they  weather  into  low  strength  clay.   The  shale  has  few  open 
Joints,  but  the  limestone  is  well  Jointed.   The  beds  of  shale  and  lime- 
stone are  essentially  horizontal. 

2.2  Shales  of  Silurian  Age 

The  Silurian  System  is  represented  in  Indiana  by  a  succession 


of  limestones  and  dolomites.   Silurian  rocks  are  exposed  at  the  surface 
in  the  southeastern  part  of  the  state.   North  of  the  Illinoian  glacial 
boundary,  glacial  drift  of  varying  thickness  covers  the  bedrock  surface. 
In  some  locations,  such  as  certain  creeks  and  river  beds,  the  glacial 
drift  has  been  removed  by  erosion  and  the  bedrock  is  exposed.   Despite 
the  predominance  of  carbonate  rocks  in  the  Silurian,  there  are  two 
formations  with  prominent  shale  lithologies,  the  "Waldron"  and  the 
"Mississinewa" . 

a)  Waldron  Shale.   This  formation  overlies  the  Laurel  Limestone 
and  is  overlain  in  turn  by  the  Louisville  Limestone.   It  ranges  from  5 
to  12  feet  in  thickness.   In  a  few  places  the  Waldron  is  entirely 
missing,  and  the  Louisville  Limestone  rests  uncomformably  on  the  Laurel 
Limestone. 

Waldron  Shale  contains  the  clay  minerals  kaolinite,  illite    and 
chlorite;  the  non-clay  fraction  consists  of  quartz,  dolomite  and  calcite. 
The  color  of  Waldron  Shale  varies  from  green  to  gray.   In  some  places  it 
is  massive  and  soft,  but  in  others,  is  fissile  and  hard.   It  is  fine 
grained  and  much  of  the  quartz,  calcite  and  dolomite  are  nearly  as  fine 
grained  as  the  clay  minerals. 

b)  Mississinewa  Shale.  This  formation  overlies  the  Louisville 
Limestone  and  is  overlain  in  turn  by  the  Liston  Creek  Limestone.  In 
places  along  the  Mississinewa  River  it  is  more  than  50  feet  thick. 

Shale  in  the  Mississinewa  is  mineralogically  similar  to  Waldron 
Shale,  except  that  the  former  commonly  lacks  kaolinite.   In  many  places 
Mississinewa  rocks  are  primarily  dolomite  rather  than  shale,  as  the 
dolomite  beds  consistently  comprise  more  than  50  percent  of  the  sequence. 


Mississinewa  Shale  is  gray  to  blue  on  fresh  surfaces  and  light 
brown  when  weathered.   The  shale  is  massive  and  bedding  is  not  apparent. 
It  is  dominantly  fine  grained,  but  in  some  places  silt  or  sand  give  the 
unit  a  silty  or  sandy  appearance. 

2. 3  Shales  of  Devonian  Age 

Similar  to  the  Silurian,  the  Devonian  System  is  also  represented 
in  Indiana  by  a  succession  of  limestones  and  dolomites.   They  are  exposed 
at  the  surface  in  southeastern  Indiana,  but  are  otherwise  covered  by 
glacial  drift  of  varying  thickness.   There  is  only  one  shale  formation, 
the  "New  Albany",  contained  in  the  Devonian  sequence. 

a)  Npw  Albany  Shale.   This  shale  is  partially  Devonian  and 
Mississipian  in  age,  since  the  unit  extends  above  and  below  the  boundary 
of  these  two  geologic  systems.   The  North  Vernon  Limestone  underlies  the 
New  Albany  Shale  and  the  Rockford  Limestone  of  Mississippian  age  overlies 
it. 

New  Albany  Shale  varies  in  thickness  from  80  to  150  feet  along 
its  outcrop.   Illite,  kaolinite  and  chlorite  are  common  clay  minerals. 
Quartz  is  the  most  abundant  non-clay  mineral  and  is  associated  with 
feldspar,  calcite,  dolomite  and  phosphate  minerals.   Pyrite  is  commonly 
present  as  coarse  to  fine  crystals  in  the  form  of  nodules  or  concretions. 

New  Albany  Shale  is  dark  gray,  dark  olive  green,  or  black,  and 
weathers  to  light  gray,  brown  or  maroon  after  a  few  years  of  exposure. 
It  is  found  in  almost  all  known  locations  as  thinly  bedded  to  fissile, 
fine  grained  shale.   Some  of  its  quartz  and  pyrite  grains  are  large 
enough  to  fall  within  the  sand  size. 


2.1*  Shales  of  Mlssissippian  Age 

Mississippian  rocks  in  Indiana  are  exposed  in  a  band  that  trends 
in  a  northwest-southeast  direction  across  the  approximate  center  of  the 
state.   The  oldest  rocks  (Kinderhook)  are  at  the  eastern  edge  of  this 
band,  and  the  youngest  rocks  (Chester)  are  at  the  western  edge.  Much 
of  the  band  of  Mississippian  rocks  is  buried  by  glacial  drift. 

Rocks  of  the  Mississippian  System  are  assigned  to  four  series  in 
Indiana:   Kinderhook,  Osage,  Meramec ,  and  Chester.   The  only  shale  of 
the  Kinderhook  Series  occurs  in  the  top  portion  of  the  New  Albany  Shale, 
which  has  been  discussed  under  the  Devonian  System.   Rocks  of  the  Osage 
Series  consist  of  shales  and  limestones.   Meramec  rocks  are  mostly 
limestones  and  dolomites  and  contain  practically  no  shale.   Rocks  of 
the  Chester  Series  are  composed  of  limestones,  sandstones  and  shales. 

a)   Osage  Shales.   Rocks  of  the  Osage  Series  are  assigned  to  the 
Borden  Group  and  are  popularly  known  as  "Borden  rocks".   The  shales  of 
this  group  occur  in  two  formations,  the  New  Providence  Shale,  which  is 
the  oldest  formation  of  the  Borden  Group  and  the  Locust  Point  Formation, 
which  lies  directly  above  the  New  Providence.   These  two  have  a  similar 
lithology  and  are  difficult  to  distinguish. 

Rocks  of  the  Borden  Group  lie  in  a  narrow  band  about  12  to  15 
miles  wide  trending  from  New  Albany,  on  the  Ohio  River,  to  Lafayette. 
From  the  Illinoian  glacial  boundary  northward,  glacial  drift  of  varying 
thickness  covers  the  Borden  Shales,  but  they  are  locally  exposed  in 
places .   The  Borden  Group  overlies  either  the  Rockford  Limestone  or  the 
New  Albany  Shale  because  of  a  prominent  unconformity  which  cuts  across 
those  geologic  units. 


Borden  3ial.es  contain  illite,  kaolinite,  chlorite,  quartz  and 
feldspar.   The  non-clay  particles  in  Borden  Shales  are  commonly  silt 
size.   The  color  ranges  from  blue  gray  to  brown.   In  most  places  these 
shales  are  massive  to  blocky  on  fresh  surfaces,  but  on  weathered  surfaces 
they  display  definite  partings  and  break  out  in  small  pieces.   The  shales 
vary  from  soft  to  very  hard. 

b)  Chester  Shales.   The  rocks  of  the  Chester  Series  consist  of 
shales,  sandstones  and  limestones.   This  series  is  more  variable  in 
mineralogy,  thickness,  and  physical  properties,  both  laterally  and 
vertically,  than  the  shale  units  previously  described. 

Chester  rocks  crop  out  in  a  band  west'  of  the  Borden  rocks.   The 
outcrop  belt  extends  from  the  Ohio  River  to  a  point  midway  between 
Indianapolis  and  Terre  Haute. 

The  Bethel  Formation  is  stratigraphically  the  lowest  formation 
of  the  Chester  Series  which  contains  shale.   The  Bethel  consists  of 
dark  gray  shales  and  argillaceous  sandstones,  5  to  30  feet  thick.   Most 
of  the  Bethel  shales  are  soft,  and  their  grain  size  ranges  from  coarse 
to  fine. 

All  shales  of  the  Chester  Series  are  variable  in  physical 
properties  and  mineralogy,  particularly  in  the  lateral  direction.   In 
many  locations  the  shale  grades  laterally  into  sandstone  or  limestone 
in  less  than  a  mile,  and  in  other  places  in  a  matter  of  tens  of  feet. 
The  dominant  clay  minerals  are  illite,  kaolinite  and  occasionally 
montmorillonite.   Quartz  is  the  main  non-clay  mineral.   Feldspar  and 
calcite  may  be  present  in  small  quantities. 
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2.5  Shales  of  Pennsylvania!!  Age 

Rocks  of  the  Pennsylvanian  System  lie  west  of  the  Mississippian 
outcrop,  in  a  belt  extending  from  the  Ohio  River  northward  to  Lafayette, 
and  then  westward  to  the  Indiana-Illinois  state  boundary.   North  of  the 
Illinoian  glacial  boundary,  glacial  drift  of  varying  thickness  covers 
most  Pennsylvanian  rocks. 

Pennsylvanian  formations  are  stratigraphically  complex  because 
of  common  changes  from  one  rock  type  to  another  over  relatively  short 
distances.   In  addition,  rocks  of  a  specific  lithologic  type  are 
similar  mineralogically  from  one  Pennsylvanian  formation  to  another, 
making  it  difficult  to  distinguish  between  the  formations  using 
lithology  alone. 

Two  types  of  shales  are  found  in  Pennsylvanian  rocks  in  Indiana: 
l)  dark-gray  to  black,  fine  grained  thinly  bedded  shale;  and  2)  light-gray 
silty  thick  bedded  shale. 

Pennsylvanian  shales  have  less  quartz  and  feldspar  than  the  shales 
previously  discussed.   The  common  clay  minerals  are  illite,  kaolinite 
and  chlorite.   They  also  contain  traces  of  iron. 

A  summary  of  the  available  data  sets  based  on  the  peolopical  descrip- 
tors is  given  in  '"able  1. 

3.   Physiography  of  Indiana 

It  was  decided  to  use  the  bedrock  physiography  as  modified  by 
Wayne  (1956)  as  the  basis  for  the  physiographic  description.  Figure 
2  shows  these  bedrock  physiographic  units. 
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TABLE  1 


DATA  SETS  IN  STORAGE  SYSTEM  ACCORDING  TO  GEOLOGICAL  PKGORTrTTON 


Geological  System 


Geological  Stage 


No.  of  Data  Sets 


Pennsylvanian 


Shelburn  Formation 
Dugger  Formation 
Petersburg  Formation 
Brazil  Formation 
Mansfield  Formation 


Total 


Mississippian 


Core  Limestone 
Palestine  Sandstone 
Watersburg  Sandstone 
Tar  Springs  Formation 
Glen  Dean  Limestone 
Hardinsburg  Formation 
Haney  Limestone 
Big  Clifty  Formation 
Elwren  Formation 
Sample  Formation 
Bethel  Formation 
Borden  Group 
Locust  Point  Formation 
New  Providence  Shale 
New  Albany  Shale 


Total 


2 

IT 
5 

1? 
3 

20 

3 

3 
1 

3 

16 
1* 
I. 

5_ 

105 


Devonian 


Antrim  Shale 
New  Albany  Shale 


Total 


1 

C_ 
7 


Ordovician 


Whitewater  Formation 
Dillsboro  Formation 
Kope  Formation 


Total 


2 
6 
_12_ 
20 
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Figure  2.   Bedrock  Physiographic  Units  of  Indiana  (after  Wayne,  1°5^) 
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In  Indiana  the  bedrock  physiographic  units  are  regional  features 
that  were  formed  by  differential  rates  of  erosion  on  gently  dipping 
Paleozoic  shales,  limestones,  and  sandstones.   These  physiographic 
units  must  have  been  in  existence  before  the  Pleistocene  epoch,  for  the 
units  are  well  defined  even  in  northern  Indiana,  where  they  were 
protected  from  extensive  dissection  during  the  Pleistocene  by  a  mantle 
of  glacial  drift  (Wayne,  1956). 

The  well-differentiated  physiographic  units  of  central  and 
northern  Indiana  which  had  been  formed  during  the  interval  that  followed 
late  Tertiary  peneplation  and  preceded  glaciation  were  covered  by  glacial 
deposits  and  preserved  from  further  modification  (Wayne,  1956). 

The  description  of  the  different  units  which  follows  is  taken  from 
Wayne  (1956). 

3.1  Dearborn  Upland 

The  Dearborn  upland,  which  lies  within  the  area  covered  by  ice 
during  Illinoian  glaciation, is  characterized  by  smooth,  steep  slopes  and 
long,  flat-topped  fingers  of  upland  between  deeply  intrenched  valleys. 
Most  of  the  valley  bottoms  are  narrow.   The  upland,  ranging  from  950  to 
1,000  feet  in  altitude,  is  a  glacially  modified  remnant  of  the  Lexington 
peneplain. 

Buried  valleys  that  have  been  traced  are  deeply  incised  below  a 
fairly  flat  upland  whose  slope  coincides  with  the  gentle  dip  of  the  rocks. 
Glacial  deposits  thicken  from  less  than  50  feet  in  Franklin,  Union,  and 
Fayette  Counties  to  more  than  200  feet  over  the  north  and  west  margins 
of  this  upland. 
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3.2  Muscatatuck  Regional  Slope 

The  back  slope  of  the  escarpment  formed  by  the  resistant  Silurian 
limestones  of  southern  Indiana  is  a  gently  westward-dipping  structural 
plain  called  the  Muscatatuck  Regional  Slope.   It  is  covered  with  thin 
Illinoian  drift  in  Indiana  and  can  be  recognized  for  a  short  distance 
north  of  the  Wisconsin  glacial  boundary. 

North  of  the  Wisconsin  boundary  the  Muscatatuck  Regional  Slope 
retains  its  characteristics  as  far  as  the  bedrock  structure  has  a 
distinct  westerly  dip.   In  Rush  and  Henry  Counties,  Middle  Silurian 
limestones  increase  in  thickness  from  about  50  feet  to  more  than  200 
feet,  and  the  regional  dip  ranges  from  about  30  feet  per  mile  to  the 
southwest  to  almost  the  vanishing  point  over  the  broad  top  of  the 
north  end  of  the  Cincinnati  Arch.   The  northern  boundary  of  the 
Muscatatuck  Regional  Slope  has  been  placed  along  the  upper  reaches  of 
a  large  bedrock  valley  in  southern  Hancock  County  and  across  Henry 
County. 

A  northward-trending  belt  of  fairly  rugged  morainic  topography 
lies  over  the  Laughery  escarpment  in  Randolph,  Henry,  and  Rush  Counties. 
The  high  bedrock  surface  along  this  escarpment  probably  contributed  to 
the  formation  of  a  re-entrant  between  ice  lobes  in  this  area.   Drift  is 
thin  and  outcrops  of  the  underlying  bedrock  are  abundant  along  the 
eastside  of  the  regional  slope,  but  glacial  deposits  increase  to  about 
200  feet  in  thickness  toward  the  west  edge  of  the  slope. 

3. 3  Scottsburg  Lowland 

The  Muscatatuck  regional  slope  grades  westward  into  a  narrow 
lowland  formed  on  Devonian  and  lower  Mississippian  shales.  Most  of  this 
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region,  the  Scottsburg  Lowland,  is  between  600  and  700  feet  in  altitude, 
and  relief  is  slight.   Thin  Illinoian  drift  blankets  this  region.   The 
Knobstone  escarpment  rises  abruptly  above  the  west  edge  of  the  lowland 
to  form  its  boundary. 

Glacial  deposits,  which  thicken  north  of  Johnson  County,  entirely 
obscure  this  lowland.   Because  the  glacial  deposits  over  this  physio- 
graphic unit  are  thick,  in  places  exceeding  U00  feet,  not  many  water 
wells  reach  bedrock.   Therefore  only  a  few  of  the  larger  drainage  lines 
have  been  identified  across  the  buried  lowland. 

Although  prior  to  the  Pleistocene  Epoch  the  Scottsburg  Lowland 
in  Indiana  was  crossed  by  five  large  rivers, none  of  these  followed  it 
for  any  great  distance.   The  Lowland  is  primarily  the  result  of  more 
rapid  erosion  on  shales  than  on  limestones  and  sandstones  which 
underlie  adjacent  uplands.  Much  of  this  Lowland  probably  was  formed  in 
late  Tertiary  and  early  Pleistocene  time,  although  degradation  continued 
in  the  southern  part  throughout  most  of  the  Pleistocene. 

The  present  topography  over  this  preglacial  Lowland  includes 
some  of  the  flattest  parts  of  the  Tipton  Till  Plain.   The  few  moraines 
are  weakly  developed,  and  postglacial  dissection  has  been  slight, 
Drift  is  so  thick  that  the  underlying  shales  are  exposed  only  along 
the  Wabash  Valley  near  Delphi,  where  postglacial  erosion  has  un- 
covered the  New  Albany  (black)  Shale. 

3.1+  Norman  Upland 

The  Knobstone  escarpment ,  by  far  the  most  prominent  topographic 
feature  of  southern  Indiana,  rises  U00  to  600  feet  above  lowland  east 
of  it  and  reaches  an  altitude  of  1,000  feet  between  the  Ohio  and 
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Muscatatuck  Rivers.   This  escarpment  forms  the  east  edge  of  the  Norman 
Upland,  a  physiographical  unit  that  is  on  the  dips]  ope  of  a  cuesta 
produced  by  resistant  sandstones  and  siltstones  of  the  upper  part  of  the 
Borden  (lower  Mississippian)  series.   The  western  boundary  of  the  Upland 
is  gradational  and  is  drawn  where  sinkholes  of  the  karst  plain  to  the 
west  become  predominant  over  valleys  that  have  been  eroded  by  running 
water. 

The  Norman  Upland  is  a  remnant  of  the  late  Tertiary  Lexington 
peneplain.   Its  altitude  exceeds  1,000  feet  in  the  south  but  decreases 
gradually  northward  to  650  ft  in  Warren  County.   The  upland  is  maturely 
dissected  in  the  unglaciated  part  of  southern  Indiana.   Dissection  of 
this  upland  also  undoubtedly  continues  northward,  but  many  of  the  smaller 
buried  ravines  probably  never  will  be  located   partly  because  ihe 
locations  of  these  buried  ravines  are  unknown  and  partly  because  the 
buried  part  of  the  Norman  Upland  has  been  less  dissected  than  that  south 
of  the  glacial  boundary. 

The  width  of  the  Norman  Upland  broadens  northwestward  from  less 
than  1  mile  at  the  Ohio  River  to  about  Uo  miles  north  of  East  White 
River  and  then  narrows  again  farther  north.   Glacial  drift  is  less  than 
50  feet  thick  over  much  of  the  Upland,  and  cut crops  are  fairly  common. 
Drift  thickens  to  more  than  200  feet , however,  over  the  intrenched  bedrock 
valleys  in  Hendricks  and  Montgomery  Counties. 

A  complex  group  of  moraines  deposited  during  a  readvance  of 
Wisconsin  ice  constitutes  much  of  the  present  surface  topography  above 
the  Norman  Upland  in  Fountain,  Montgomery,  and  Hendricks  Counties. 
Although  the  buried  Knobstone  escarpment  has  no  recognizable  surface 
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expression  so  far  north,  the  higher  bedrock  must  have  influenced  the 
deposition  of  these  moraines. 

3.5  Mitchell  Plain 

The  Mitchell  Plain  is  a  karst  terrain  that  was  formed  upon  the 
middle  Mississippian  limestones  in  southern  Indiana.   A  few  miles  north 
of  the  Wisconsin  boundary  these  limestones  are  covered  by  an  overlap  of 
lower  Pennsylvania^  elastics.   A  few  karst  features  can  be  seen  through 
thin  drift  in  Putnam  County  and  the  buried  Mitchell  Plain  exists  as  far 
north  as  southwestern  Montgomery  County,  where  solution-widened  bedding 
planes  and  fractures  are  reported  occasionally  by  well  drillers.   No 
other  karst  features,  however,  have  been  recognized  beneath  the  drift, 
which  ranges  from  the  vanishing  point  to  150  feet  in  thickness. 

3.6  Crawford  Upland 

The  most  rugged  topography  of  southern  Indiana  is  found  in  the 
Crawford  Upland,  a  physiographic  unit  that  is  separated  from  the  Mitchell 
Plain  to  the  east  by  the  rugged  Chester  escarpment.   Trenchlike,  flat- 
bottomed  valleys,  rock  benches,  and  local  structural  plains  typify  its 
land  forms.   Caverns  and  karst  valleys  are  common,  and  remnants  of  the 
Lexington  peneplain  are  recognizable  at  some  places.   The  western 
boundary  of  the  Crawford  Upland  is  located  where  the  angular  topography 
formed  on  alternating  massive  sandstones,  shales,  and  limestones  of  the 
Chester  (upper  Mississippian)  series  grade  into  rounded  land  forms  that 
are  a  result  of  erosion  on  Pennsylvanian  shales  and  coals. 

North  of  the  Wisconsion  drift  margin  the  upland  widens  somewhat 
to  include  the  region  of  gorgelike  valleys  cut  into  the  Mansfield  sand- 
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stone  in  Parke,  Vermillion,  Fountain,  and  Warren  Counties.   Relief 
features  characteristic  of  the  Crawford  upland  continue  as  far  north 
as  the  bluff  overlooking  the  Teays  Valley.   Part  of  the  preglacial 
divide  between  the  Teays  and  Wabash  drainage  basins  was  in  northern 
Parke  County, and  Sugar  Creek  now  is  trenching  part  of  this  former 
upland. 

During  part  of  the  Wisconsin  glaciation,  the  north  end  of  the 
Crawford  Upland  lay  between  ice  lobes  in  central  Indiana  and  east- 
central  Illinois.   Although  it  is  the  location  of  many  rugged  morainic 
features,  drift  is  thin  over  most  of  this  preglacial  upland,  and  many 
present  streams  have  cut  valleys  in  bedrock.   Thick  drift  is  limited 
largely  to  buried  valleys. 

3.7  Sullivan  Lowland 


Prior  to  glaciation  a  topography  of  subdued  land  forms  which 
had  moderate  relief  along  the  major  valleys  was  formed  upon  the 
Pennsylvanian  strata  of  central  Illinois  and  southern  Indiana,  where 
weak  shales  permitted  valley  widening  to  a  greater  extent  than  was 
possible  in  areas  of  more  resistant  rock.   Most  of  the  bedrock  slopes 
are  gentle,  smooth  inclines,  unlike  the  angular  topography  of  the 
Crawford  Upland,  and,  except  near  the  bedrock  strata  of  major  valleys 
such  as  the  Wabash,  deep  stream  courses  are  uncommon.   Because  the 
bedrock  surface  in  Sullivan  County,  Ind.,  seems  to  exemplify  the 
characteristics  of  this  unit,  the  name  "Sullivan  Lowland"  was  used  for 
it  by  Wayne  (1956). 

Little  of  the  Sullivan  Lowland  lies  north  of  the  Wisconsin  boundary 
in  Indiana.   It  is  restricted  to  southern  Parke  County,  Vermillion  County, 
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and  the  west  edge  of  Warren  County,  where  it  seems  to  merge  with  the 
Knobstone  escarpment  along  the  south  bluff  of  the  Teays  Valley  near  the 
state  line. 

3.8   Bluffton  Plain 

A  nearly  flat  limestone  upland  that  slopes  gently  to  the  north 
constitutes  the  bedrock  surface  in  east-central  Indiana.   It  is  bounded 
by  a  large  preglacial  valley  from  northern  Allen  County  to  near 
Logansport,  in  Cass  County.   On  the  west  it  grades  into  the  Scottsburg 
Lowland  and  on  the  south  lies  next  to  the  Muscatatuck  Regional 
Slope  and  the  Dearborn  Upland.   The  boundary  between  this  upland  plain 
and  the  Muscatatuck  Regional  Slope  is  a  transitional  one  in  Hancock 
and  Henry  Counties.   It  has  been  drawn  across  a  narrow  divide  between 
the  Whitewater  Bedrock  and  Wabash  Bedrock  drainage  basin  in  south- 
western Henry  County  and  then  follows  the  course  of  a  bedrock  valley 
westward  across  Hancock  County.   Because  the  flat,  plainlike  character 
of  this  upland  is  well  developed  in  Wells  County  and  drift  is  thin 
enough  near  Bluffton  to  expose  the  underlying  dolomites,  the 
name  "Bluffton  Plain"  was  proposed  for  this  region. 

This  plain  is  about  900  feet  above  sea  level  in  Randolph  and  Jay 
Counties.   Its  altitude  decreases  to  750  feet  near  Teays  Valley  in 
Blackford  and  Grant  Counties  and  to  650  to  TOO  feet  in  Wabash  County. 
North  of  Teays  Valley  it  never  rises  much  above  750  feet  in  altitude 
but  slopes  into  a  lowland  to  the  north.   The  plain  was  formed  upon 
the  thick  sequence  of  Silurian  limestones  and  dolomites  of  northern 
Indiana,  and  its  slope  corresponds  closely  to  the  regional  dip  on  the 
north  end  of  the  Cincinnati  Arch. 
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Karst  phenomena  were  at  least  moderately  well  formed  upon  the 
limestone  plain  prior  to  its  burial.   Discoveries  of  caverns  have  been 
reported  by  water  well  drillers  in  Adams,  Wells,  and  adjacent  counties 
near  Teays  Valley  and  some  of  its  deeply  incised  tributaries.  Some  of 
these  cavernous  openings  in  the  dolomite  are  partly  filled  with  red  silt 
and  clay;  others  are  open  and  yield  large  quantities  of  ground  water. 
Small  dolines  may  be  observed  in  some  of  the  reefs  exposed  near  Lagro , 
in  Wabash  County. 

3.9  Rensselaer  Plateau 


The  northeastward-sloping  upland  of  northern  Illinois  seems  to 
extend  into  Indiana.   It  was  named  the  "Rensselaer  plateau,"  because 
the  features  of  this  unit  are  well  displayed  in  the  bedrock  topography 
of  Jasper  County,  Ind.,  in  the  vicinity  of  Rensselaer.   This  unit  is 
largely  a  dip  slope  that  was  formed  upon  Middle  Silurian  dolomites 
and  limestones,  but  Devonian  and  Mississippian  limestones  also  underlie 
it. 

To  the  south,  the  Rensselaer  plateau  ends  at  Teays  Valley,  and  its 
northeastern  boundary  is  placed  along  a  large  valley  in  LaPorte  and 
Marshall  Counties. 

Although  glacial  deposits  are  less  than  50  feet  thick  over  a  large 
area  in  the  Kankakee  Valley,  no  bedrock  outcrops  are  known.   Drift 
thickens  to  nearly  300  feet  in  the  large  valley  at  the  northeastern 
boundary  of  the  region.  Much  of  the  bedrock  surface  has  so  little 
relief  that  thickening  of  drift  along  the  Valparaiso  moraine  showed 
up  readily. 
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3.10  DeKalb  Lowland 

The  bedrock  physiography  of  most  of  the  northern  two  tiers  of 
counties  in  Indiana  is  a  broad  lowland  formed  upon  Upper  Devonian  and 
Lower  Mississippian  shales.  Altitudes  range  from  ^00  feet  in  valleys 
beneath  the  drift  in  DeKalb  and  LaPorte  Counties  to  nearly  TOO  feet  in 
St.  Joseph  and  Elkhart  Counties.   Rocks  underlying  this  lowland  dip 
gently  northward  toward  the  Michigan  Basin.  The  term  "DeKalb  Lowland" 
was  proposed  for  this  physiographic  upland.   It  was  named  from  DeKalb 
County,  where  the  Lowland  is  well  formed  and  where  some  of  the  thickest 
glacial  deposits  in  Indiana  are  found. 

Because  of  the  great  thickness  of  drift  and  the  abundance  of 
buried  gravel  beds,  few  wells  drilled  over  this  Lowland  reach  bedrock. 
A  notable  exception  is  along  the  St.  Joseph  Valley,  where  several  test 
holes  in  the  industrial  areas  of  South  Bend,  Mishawaka,  and  Elkhart  have 
reached  bedrock. 

The  topography  of  the  DeKalb  Lowland  should  have  only  moderate 
relief  except  along  intrenched  valleys  like  the  topography  of  the 
Scottsburg  Lowland  of  central  and  southern  Indiana,  which  is  formed 
upon  rock  similar  in  lithology  and  structure. 

The  minimum  and  maximum  thicknesses  of  glacial  drift  recorded 
in  this  region,   namely,  28  feet  in  Mishawaka  and  525  feet  just  west 
of  Elkhart,  are  within  a  few  miles  of  each  other.  Drift  less  than 
100  feet  thick  is  exceptional,  and  most  of  the  region  is  covered  with  200 
to  350  feet  of  Pleistocene  deposits.   Some  of  the  greatest  local  surface 
relief  on  glacial  deposits  in  Indiana  is  over  the  DeKalb  Lowland. 
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3.11  Lake  George  Upland 

The  upper  part  of  the  lower  Mississippian  shales  and  sandstones 
at  the  south  margin  of  the  Michigan  Basin  produced  an  abrupt  escarp- 
ment which  crosses  the  northeast  corner  of  Indiana.   The  Coldwater 
Formation  of  early  Mississippian  age  seems  to  contain  the  resistant  beds 
that  form  the  crest  of  this  scarp.   The  name  "Lake  George  Upland"  was 
proposed  for  the  physiographic  unit  bounded  on  the  south  by  the  escarp- 
ment and  extending  northward  on  the  dip  slope.   Lake  George  in  Steuben 
County,  Ind.,  and  Branch  County,  Mich.,  is  within  a  group  of  moraines 
that  overlies  this  upland.   The  bedrock  surface  beneath  Lake  George  is 
750  to  800  feet  in  altitude  and  is  buried  beneath  200  to  250  feet  of 
drift. 

Altitudes  in  the  northern  part  of  Steuben  County  show  that  the 
upland  is  between  750  and  900  feet  above  sea  level  and  stands  as  much 
as  U00  feet  above  the  DeKalb  Lowland  south  of  it.   Glacial  drift  thins 
abruptly  over  the  scarp  from  more  than  500  feet  to  150  to  200  feet. 
Surface  altitudes  exceed  1,200  feet  in  this  part  of  Indiana,  and  some 
of  the  steepest  morainic  topography  of  the  state  is  found  here. 
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4.  Tests  Performed  on  Shale  Samples. 

A  number  of  tests  are  performed  on  shale  samples  submitted  to  the 
laboratory.   These  tests  must  comply  with  standard  testing  methods  in 
order  to  use  the  data'  in  any  statistical  analysis.   The  main  contributor 
of  data  to  the  present  system  is  the  Indiana  State  Highway  Commission. 
Results  of  their  testing  are  reported  on  a  standard  form,  a  reduced 
copy  of  which  is  shown  in  Figure  3.   The  different  testing  methods  used 
by  ISHC  will  be  discussed  in  this  section. 

The  General  Physical  Description  (see  Figure  3)  of  the  material 
is  based  on  the  reporter's  subjective  opinion  and  will  not  be  discussed 
any  further. 

A  number  of  the  tests  are  performed  according  to  AASHTOstandard 
procedures.   The  designator  numbers  of  these  tests  are  given  in  Table  2. 
The  tests  not  covered  in  this  table  are:  pH,  loss  of  ignition,  soil 
textural  classification  and  all  the  tests  for  shale  classification. 
The  soil  pH  is  determined  by  means  of  an  electrical  pH  -  meter  in  a  1:1 
soil-water  suspension.   The  method  is  described  in  Part  II  of  ASTM 
D  2976  (1972).   The  textural  classification  is  done  according  to  Figure  4, 

Tests  to  determine  the  loss-on-ignition,  slake  durability  index, 
slaking  index,  fissility  number,  modified  soundness  and  classification 
of  shale  are  described  in  Appendix  A. 

5.  Storage  of  Data. 

Laboratory  test  data  on  shales  is  reported  by  the  Indiana  State 
Highway  Commission  on  a  standard  form,  a  copy  of  which  is  shown  in 
Figure  3.   For  the  purpose  of  this  study  it  was  decided  not  to  include 
all  the  data  on  Figure  3.   The  data  of  importance  to  this  study  are  the 
description  of  the  sample  and  the  laboratory  testing  results.   It  was 
therefore  decided  to  include  only  those  data  indicated  by  a  bar  on  Figure 
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Table  2 
AASHTO  Designator  Numbers  for  Standard  Tests  Used  by  ISHC 

Title  of  Test  Designation 

1.  Dry  Preparation  of  Disturbed  Soil  and  T87-72 

Soil  Aggregate  Samples  for  Test 

2.  Wet  Preparation  of  Disturbed  Soil  Samples  T186-49 

for  Test 

3.  Particle  Size  Analysis  of  Soils  T88-72 

4.  Determining  the  Liquid  Limit  of  Soils  T89-68 

5.  Determining  the  Plastic  Limit  and  Plasticity  T90-70 

Index  of  Soils 

6.  Determining  the  Shrinkage  Factors  of  Soils  T92-68 

7.  Moisture-Density  Relations  of  Soils  Using  a  T99-74 

5.5  lb.  Rammer  and  a  12-in.  Drop 
(methods  B  &  D) 

8.  Specific  Gravity  of  Soils  T100-74 

9.  The  California  Bearing  Ratio  T193-72 

10.  Density  of  Soil-in-Place  by  Block,  Chunk  or  T233-70 

Core  Sampling 

11.  Classification  of  Soil  and  Soil-Aggregate  M145-73 

Mixtures  for  Highway  Construction  Purposes 
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Figure  h.      Details  on  Determining  the  Textural  Classification  of 
Shale  Samples 
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5.   The  original  data  sheets  should  be  consulted  if  any  other  data  are 
required. 

The  data  had  to  be  coded  for  computer  purposes.   A  four  card  series 
is  used  to  put  all  the  data  of  Figure  5  into  storage.   The  format  of 
the  coding  is  given  in  Table  3. 

In  order  to  assess  which  data  points  are  available  in  each  data 
set,  Table  4  was  compiled.   The  crosses  indicate  the  data  points 
available.   The  descriptor  used  in  this  table  to  identify  the  samples 
is  the  laboratory  number,  which  is  given  here  in  ascending  order.   The 
data  broken  down  into  the  different  geological  units  are  given  in 
Appendix  B. 

Because  of  the  small  number  (163)  of  data  sets  presently  available, 
a  number  of  parameters  are  given  as  words,  e.g.,  all  the  data  on  card  1 
and  the  two  classifications  on  card  3  (see  Table  3).   The  alternative 
was  to  use  numbers  as  codes  for  these  parameters  with  'word-files',  in 
order  to  get  a  print-out  with  the  correct  descriptions  on  it.   This 
alternative  is  not  justified  at  this  stage.   The  data  sets  can,  however, 
be  changed  easily  if  this  feature  is  necessary  in  the  future. 

The  data  sets  can  also  be  enlarged  in  the  future,  i.e.,  more  than 
four  cards  can  be  used  for  each  data  set.   This  will  only  be  justified 
if  more  data  on  other  parameters  became  available  for  a  large  number  of 
data  sets. 

The  data  stored  in  the  present  system  are  based  on  that  provided 
by  Indiana  State  Highway  Commission  as  well  as  the  results  of  studies 
on  shales  at  Purdue  University.   Table  5  gives  the  laboratory  numbers 
of  the  samples  used  by  Deo  (1972),  Chapman  (1975),  and  Bailey  (1976). 
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Figure  5.   Identification  of  Data  Included  in  This  Study 
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Table  3 


Coding  of  Four  Card  Series 


Parameter 


Symbols 
(for  program  and 
tables  and  figures) 


Format 
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Card  1 

Laboratory  number 

NO 

110 

1  - 

10 

County 

CO 
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11  - 

20 

Geological  system 

GSY1.GSY2 

A10A5 

21  - 

35 

Geological  series 

GSE1,GSE2 

A10A5 

36  - 

50 
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F5.1 

31  - 

35 

200  rev.  soaked 
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26  - 
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QL 

F5.1 
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^5 
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F5.1 

46  - 
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Percent  sand 
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15 

51  - 

55 
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15 

56  - 

60 

Percent  clay 
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15 

61  - 
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15 

66  - 

70 

30 


Table  3  (Continued) 


Parameter 


Symbols 
(for  program  and 
tables  and  figures) 
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Laboratory  number  NO 
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Natural  moisture  content  AMC 
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Moisture  density*  MD 
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Maximum  wet  density  WD1 
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As  Compacted  CBR  CBR1 

After  Soaking  CBR  CBR2 

Average  %  Swell  AS 
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F5.1 

52  - 

55 

F5.1 
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*Coding  for  moisture  density: 

1  -  Test  done  on  minus  No.  4  material 

2  -  Test  done  on  minus  3/4  inch  material 

3  -  When  there  is  a  special  note  on  the  laboratory  report  which 

should  be  referred  to 

4  -  Nothing  indicated  on  test  report 
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TABLE  5 


Laboratory  Numbers  of  Shale  Samples  Utilized  by 
Researchers  at  Purdue 


Deo  (1972) 
Samples  1  to  15 

Chapman  (1975) 
Number  in  report 
1 

2 
3 
4 
5 
6 


ISHC  Laboratory  No, 
73-51703 
74-54621 
74-54684 
74-54716 
75-55018 
75-55315 


Bailey  (1976) 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 


74-54878 
74-54684 
75-55653 
74-54716 
74-54836 
75-55044 
74-54767 
73-51703 
74-54973 
74-54972 
75-55564 
75-55316 
75-55315 
75-55731 
75-55505 
74-54621 
75-55718 
75-55486 
75-55487 
75-55018 
75-55291 
76-55014 
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There  are  some  discrepancies  between  the  values  of  the  parameters 
as  given  in  the  report  by  Chapman  (1975)  and  the  values  given  by  ISHC 
on  their  results  sheet.   The  different  values  are  given  in  Table  6. 
In  order  to  limit  the  variables  contributing  to  differences  in  test 
parameter  values  it  is  necessary  to  use  data  from  as  few  different 
sources  as  possible.   With  this  in  mind  it  was  decided  to  use  the  values 
supplied  by  ISHC. 

6.   Retrieval  of  Data 

There  are  basically  three  purposes  for  the  retrieval  of  data: 

a)  To  get  print-outs  of  all  the  data  available. 

b)  To  get  print-outs  based  on  some  descriptor  in  some 
order, e.g.,  geological  stage  in  alphabetical  order. 
The  printed  data  sets  belonging  to  a  certain 
geological  stage  can  then  be  separated  for  study. 

c)  To  use  certain  sets  of  data  for  statistical  manipulation. 
The  program  for  printing  out  all  the  data  sets  in  the  order  which 

they  are  on  the  data  cards  is  given  in  Appendix  C.   This  program  was 
written  for  the  CDC  system  of  the  Purdue  University  Computer  Center 
and  might  need  changes  in  order  to  use  it  on  any  other  system.   A 
typical  print-out  of  a  data  set  is  given  in  Figure  6. 

The  sorting  of  the  data  is  accomplished  by  using  the  SORT/MERGE 
package  on  the  CDC  system.   Sorting  can  be  done  in  an  ascending  or 
descending  order  based  on  any  parameter.   Before  the  sorting  can  be 
done,  each  data  set  must  be  changed  to  a  string  of  320  characters  in  the 
computer,  and  after  sorting  it  is  broken  up  again  into  4  groups  of  80 
characters  each.   Printing  can  now  take  place  by  using  the  program 
given  in  Appendix  C.   Detail  of  this  is  given  in  Table  7.   Detail  of 
the  sorting  is  also  given  in  this  table. 


w 
-J 

2 


u 
x 

CO 


T3 
01 
T3 
•H 
> 
O 

CM 

QJ 

in 

o 

-o 

a 


c 

M 

0 

ex 

O 
J3 


U 

o 

QJ 


c 

0) 

> 

■H 

e> 
u) 

0) 
3 

.-I 
rt 

> 

C 

QJ 
QJ 

2 

j-j 

QJ 

M 

on 

0) 
•H 

a 

c 

rt 

a 
a) 

ij 

CJ 
CO 


U 
33 
CO 


O     (U 

■h  -a 

CO    M 

rt 
H 


CJ 

CO 


•H  -H 

3  0 

•H  _. 

►J 


u 

EC 

w 


4J  -H 

CO  0 

n)  -H 

rH  r-1 


■a 

QJ  CJ 

^i  X 

rt  en 

o  H 
to 

o  c 

o  rt 

x  m  0 

QJ  O. 

■a  rt 

C  J2 

H  CJ 


•H  CJ> 

H  X 

•H     >»  CO 

B    (J  H 

rt  a 
u 

3  O  CI 

a  o  « 

u->  0 

cu  a 

«s  XI 

rH  CJ 
CO 

■a 

Q)  CJ 

^!  3C 

rt  co 

O  H 

CO 


a 
rt 
1 
a. 
rt 
J3 
CJ 


CJ 

35 
>,  CO 
M  M 
Q 


<N     0 

a 
rt 
.r: 

o 


u 
o 

•U  IJ 

rt  cu 

u  xi 

°  § 

£>  3 

rt  Z 
►J 


Ul 


o 

-f 

CI 

• 

■ 

* 

\o 

^D 

m 

r-t 

iH 

rH 

OJ 

B 

d 

r-~ 

'-  l 

<r 

c/i 

in 

rH 

CO 

u-i 


CO 


in 
o 


m 


O 


en 


ci 


o 


CO 

O 


oj 


QJ 

g 
CO 


en 

CTi 


ON 


CO 


CO 
ON 


en 

CTn 


m 


o 


I-. 


en 
en 


in 

VD 


CO 


00 

en 


cn 
cn 


-3- 

m 


CO 

m 


oo 

OJ 


OJ 

CI 


CO 
Ov 

en 


oo 
en 


en 
en 


o 

no 


oo 
en 


l£J 


CO 
en 


CO 

en 


en 
en 


o 

en 


I** 

en 


o 
oo 


CO 

en 


rH 

<r 

sO 

rH 

m 

OJ 

oo 

rH 

O 

rH 

^O 

vO 

r^ 

m 

ro 

<r 

<r 

<r 

m 

m 

m 

m 

m 

rH 

m 

■* 

<r 

<r 

m 

m 

r- 

o- 

r» 

r-» 

o~ 

k2 


LAPORAIQRY    ML>fiHF.]<    -    -    -    -.  r  _-..-..------..-. _74Ji(*Afa£. 
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GEOLOGICAL     SYSTEM HlSSISSU't'lAN 

GCOLuCICAL    :>IRIF        _--.-------    CHESTFR 

glolocical  s'AGt    -  r  r.r„r._r„t..t— "—".-"— kLpjiC. 


"  PM  y¥  lo'C  R  A  PhlV    UN  11       ----------    CRAWFORD    Of  lam- 


slaking    INL-t'X    CYCLE    1---------         ?.£, 

S  L  A  K 1Mb     INDEX     CYCLt    S---------       51.1 


SLAKI'    Dul-AblLUY    XiNjOLX    £00    KtV    UKY       -    -       71.5 

SLAKt    DLUvAblL  11  y     iMLtX    £00    RLV    SUAKED_-_ 1 2_._b_ 

"¥CXKLr"tVrrAbTLTfy"l"SL"LX~B"du  "htV    LjHY       -    -       73.  *• 
SLAKE     U1JF-AG1L  1TY    lJ^ktX__5jJ0._«ty._EjaARE_D__=__i»7_.Ji. 

Ti¥si1Ci^Y~«U'''KtB"~~"~-"~  --------        -o 

hODIFIED    SQijMjNESS    TEST     (PtKL    LOSS)     -    -    1CP.0 


1EXTUKAL    CLASSIFICATION    -_-__-   .t„--_1— "— Jl AliJ  X_CLAY_ 

AA"SHfo""CLASsTf  I~C~ATTgn~~-~~-  "-    -""-"-    -    -    -  A    6  Ufa  I 

PLASTIC    LIMl     (I'ERC) -_  --'.— Z-"    ".." .  "    "       g*'9 

"Li"jjiti~L"l^Tf  "jPthCl    ----------       36.6 

PLASTICITY     IHiEjt     (PLRC)     --------        iq.7 


SAUL)     (PFT.C)     --------------  A 

s i l t  i£f.-i'-?J— r„r..r„r..r._r._t— "— t— t_i— "_ r..r. A*. 

"CLAY     (PF.mCT   --------------  ?3 

colloils  _ih.6iLCj„r--r-.:__t---— _— ---".: "."—"—  ~ 2JL 


NATURAL  '.jET  LENSlTY  (LH/CUFI)  -----  1M.H 


NATUPAL     CRY     LLNSlTY     <LH/CU>'>     -----     Ul,l 
N  A  Tl'HAL_  NO  lST_Lk_F_CO]^T_LN_T._(_t:t._KCJ.__-  __-„-„-_ 7_._0_ 

"SpXCTfTc   GRAVITY      - "- ~~    --------      ?,73 
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5HRINKAGE  LI^IT  (PE.RO   --------   jfc.b 

LIUf.AL    SHRIf^AGE     (Ptkl)     --------  q.7 
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Figure  6.      Typical  Print-out  of  a  Data  Set 
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TABLE  7 


Detail  of  Cards  for  Sorting  Data  on  CDC  of 
Purdue  University  Computer  Center  by 
Utilizing  the  SORT/MERGE  System 


Cards 

JOB  CARD 

PASS  = 

FORTRAN. 

L0AD(LG0,RUNL1B2) 

EXECUTE (.INPUT) 

REWIND (TAPE2) 

RETURN, LGO. 

SORTMRG. 

REWIND (TAPE3) 

FORTRAN . 

LOAD (LGO,  RUNL1B2) 

EXECUTE ( , TAPE3 , FILE4 ) 

RETURN, LGO. 

REWIND (FILE4) 

MNF(N) 

LGO(FTLE4) 

7/8/9 

PROGRAM  IN(TAPE1,TAPE2) 

INTEGER  ARRAY (320) 
1    CALL  READS(1,ARRAY(1),80,K) 

IF  (K.LE.O)  STOP 

CALL  READS  (1,ARRAY(81)80,K) 

IF  (K.LE.O)  STOP 

CALL  READS(1,ARRAY(161),80,K) 

IF  (K.LE.O)  STOP 

CALL  READS (1, ARRAY (241), 80, K) 

IF  (K.LE.O)  STOP 

CALL  WRITES (2, ARRAY, 320) 

GO  TO  1 

END 
7/8/9 

Data 
7/8/9 

SORTQ, 1,330) 
FILECTAPE2,S,D) 
KEY (A, C, 251, 5) 
FILE (TAP E3,0,D) 
RECORD CI, U, 330) 
END 


Comments 


Control.  Start  in  Column  1, 


Program  to  write  4  card 
series  (80  characters  each) 
as  one  string  of  320 
characters.   The  use  of  this 
special  input/output 
capabilities  is  described  in 
PUCC  publication  IO-XFORTIO 
Punching  starts  in  Column  7. 


Sorting  of  data,  see  description 
below.   Punching  starts  in 
Column  1.   See  PUCC  publication 
M2- SORTMRG. 
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TABLE  7  (Continued) 


7/8/9 

PROGRAM  OUT(TAPEl,TAPE2) 

INTEGER  ARRAY (320) 
1    CALL  READS (1, ARRAY, 320, K) 

IF  (K.LE.O)  STOP 

CALL  WRITES(2,ARRAY(1),80) 

CALL  WRITES(2,ARRAY(81),80) 

CALL  WRITES (2, ARRAY (161), 80) 

CALL  WRITES(2,ARRAY(241),80) 

GO  TO  1 

END 
7/8/9 

Program  of  Appendix  A 
6/7/8/9 

Sorting  of  data 

The  only  card  to  be  changed  here  is  the  KEY-card 


Program  to  write  320 
character  string  into 
4  groups  of  80  characters 
each.   See  publication 
above  for  detail.   Punching 
starts  in  Column  7. 


One  key  card  must  be  specified  for  each  field  to  be  sorted.   The  order 
of  importance  of  the  keys  is  indicated  by  the  order  in  which  the  KEY 
control  cards  occur  in  the  SORT/MERGE  control  deck. 


KEY(pl, 


,p4) 


Pi 
P2 


P3 


P4 


Sort  order 

A  =  Ascending 

D  =  Descending 

Key  mode  designator 

C  =  character  mode  key 

X  =  Fixed  point  key 

F  =  Floating  point  key 

L  =  Logical  key 

N  =  Numeric  BCD  Field. 

Position  indicator 

For  C  and  N  key  modes;  position  of  the  first  6-bit  byte  of  the  key 
L  key  mode;  relative  bit  position  of  the  first  bit  of  the  key 
X  and  F  key  modes;  relative  word  position  of  the  key 
C  and  N  key  modes;  length  of  key  in  6-bit  bytes 
L  key  mode;  length  of  key  in  bits 
X  and  F  keys;  ignored. 


For 
For 
For 
For 
For 


The  example  included  above  therefore  means  that  the  data  sets  will  be 
sorted  in  an  ascending  order  (pl=A) ;  we  are  dealing  with  character  mode 
key  (p2=C);  the  sorting  will  start  in  Column  251(p3=251)  therefore  in 
Column  11  on  card  4  (see  Table  3)  that  is  natural  wet  density;  the  sorting 
will  be  done  for  the  next  5  columns  (p4=5).   This  example  will,  therefore 
sort  the  data  according  to  values  of  natural  wet  density  in  ascending  order 
before  printing  the  data  sets  as  in  Pig.  6  . 
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Any  parameter  can  be  used  for  sorting  as  Indicated  in  Table  7. 
The  print-out  is  still  the  same  as  that  in  Figure  6  and  all  the  data 
sets  included  in  the  data  submitted  with  the  program  will  be  printed. 
Separate  parameters  can  therefore  be  sorted  but  no  detailed  sorting,  e.g.  , 
all  data  sets  containing  liquid  limits  higher  than  25%  and  belonging  to 
the  Mississippian  system  included  in  the  Clore  stage,  can  be  done  with 
the  SORT-MERGE  package.   It  was  felt  that  a  general  sorting  program  is 
not  justified  at  this  stage  because  of  the  small  number  of  data  sets  in 
the  present  system. 

A  general  sorting  program  would  be  necessary  to  use  the  stored 
data  for  statistical  manipulation.  For  the  present  study,  detailed 
sorting  was  done  by  hand  prior  to  some  of  the  statistical  analyses. 

The  program  in  Appendix  C  is  a  very  basic  program  and  is  only  meant 
for  reading  and  printing  data.   Modifications  of  this  program  would  be 
necessary  if  anything  more  is  necessary. 

7.   Statistical  Analysis  of  Data. 

The  statistical  analysis  was  done  by  using  the  routines  available 
in  the  statistical  package  for  the  social  sciences  (SPSS)  (Nie, 
et  al.,  1975).   This  package  is  available  at  the  Purdue  University 
Computer  Center. 

This  section  will  describe  the  different  analyses  performed  as 
well  as  the  results  obtained. 

7.1  Frequency  Analysis 

The  first  step  in  the  statistical  analysis  was  to  do  a  complete 
frequency  analysis  of  all  the  data  available.   For  this  analysis  the 
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CODEBOOK  routine  of  SPSS  was  used.   The  data  can  be  analyzed  as  'single 
points'  or  as  grouped  data.   The  data  were  first  analyzed  as  single 
points  in  order  to  obtain  minima  and  maxima  of  each  parameter.   The  values 
of  each  parameter  were  then  divided  into  10  intervals  in  order  to 
obtain  histograms.   A  number  of  statistics  are  computed  by  this  program. 
The  different  statistics  and  the  way  that  they  are  calculated  in  the 
program  are  given  in  Table  8. 

The  results  of  the  statistics  of  the  different  parameters  are  given 
in  Table  9.   The  first  line  of  results  for  every  parameter  gives  the 
values  obtained  by  using  single  points,  while  the  second  line  gives  the 
results  for  grouped  data.   The  histograms  are  given  in  Figure  7. 

The  last  two  columns  of  Table  9  give  the  number  of  values  which 
are  valid  and  also  the  missing  values,  e.g.!  for  SDI1,  131  data  sets  had 
values  reported  for  this  parameter  while  it  was  not  reported  for  32 
data  sets.   The  lowest  number  of  valid  values  is  26  for  the  modified 
soundness  test,  a  reasonably  low  number  of  test  values  is  also  available 
for  the  California  Bearing  Ratio  and  average  percentage  swell  (53). 
The  number  of  valid  values  will  play  a  very  important  role  when  the 
regression  analysis  is  done. 

The  low  values  obtained  for  the  skewness  for  liquid  limit  (QL) 
and  percent  silt  (PSI)  indicate  that  the  distribution  of  these  two 
parameters  are  almost  symmetric-bell  shaped  (refer  to  Figure  7). 
Low  values  of  coefficient  of  variance  were  also  obtained  for  these  two 
variables  which  confirms  the  small  relative  spread  of  values.   Low 
values  of  the  coefficient  of  variation  are  to  be  expected  for  SG  and  pH 
while  the  low  values  for  the  coefficient  of  variation  obtained  for  the 
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TABLE  8 
Statistics  Available  from  CODEBOOK 

1.  Minimum  value 

2.  Maximum  value 

3 .  Range 

4.  Mode  -  value  which  occurs  most  often 

5.  Median  -  the  numerical  value  of  the  middle  case  or  the  case  lying 
exactly  on  the  50th  percentile,  once  all  the  cases  have  been 
ordered  from  highest  to  lowest.   In  computing  the  median  from 
grouped  data,  the  values  within  each  category  are  assumed  to  be 
distributed  evenly  throughout. 

6 .  Mean 


X  = 


i=l 


X. 
l 


Variance 


N 


s2  =  i=1 


y  x.  -  n  x 

,   1 


N 


8.  Standard  deviation  -  S 

9.  Standard  Error:  If  we  were  to  draw  an  infinite  number  of  equal- 
sized  samples  from  a  given  population,  the  mean  of  each  sample 
would  be  an  estimate  of  the  true  population  mean,  but  not  all  of  them 
would  be  identical.   The  pattern  of  these  means  would  actually 
constitute  a  normal  distribution  and  would  have  a  standard  deviation. 
The  standard  deviation  of  this  distribution  is  the  standard  error. 
Standard  error  cannot  be  computed  exactly,  but  it  can  be  estimated 

by  dividing  the  standard  deviation  by  the  square  root  of  the  number 
of  cases. 
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TABLE  8  (Continued) 


10. 


Skewness 


|Li=l       li=l 


-2 
+   3X 


N 
i-1   X 


]/■} 


-x3 


{[ 


N 

i=i i 


-NX2 


J/N-l} 


3/2 


Skewness  is  a  statistic  to  determine  the  degree  to  which  a 
distribution  of  cases  approximates  a  normal  curve,  since  it  measures 
deviations  from  symmetry.   It  will  take  on  a  value  of  zero  when  the 
distribution  is  a  completely  symmetric  bell-shaped  curve.   A 
positive  value  indicates  that  the  cases  are  clustered  more  to  the 
left  of  the  mean  with  most  of  the  extreme  values  to  the  right.   A 
negative  value  indicates  clustering  to  the  right. 


11.      Kurtosis 


.in 


4X 


N 

Ix 

i=l 


+  6X2 


N 

I   X 
i=l    X 


4X~ 


[  14  <]  •  '■ 


If 


flx2 

1=1 


-NX2 


J7H 


Kurtosis  is  a  measure  of  the  relative  peakedness  or  flatness  of 
the  curve  defined  by  the  distribution  of  cases.   A  normal  distribution 
will  have  a  kurtosis  of  zero.   If  the  kurtosis  is  positive,  then  the 
distribution  is  more  peaked  (narrow)  than  would  be  true  for  a  normal 
distribution,  while  a  negative  value  means  that  it  is  flatter. 

12.   Coefficient  of  Variation 


V  =  100 


X 


This  is  a  measure  of  dispersion. 
SPSS  but  was  calculated  by  hand. 


This  value  is  not  given  by 
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different  densities  reflect  the  high  average  values  and  small  spread 
of  these  parameters.   The  relative  low  value  of  coefficient  of 
variation  for  percent  silt  (25%)  may  be  the  result  of  the  preparation 
for  the  grain-size  analysis.   The  method  used  can  result  in  breaking  down 
the  material  to  a  reasonably  constant  size  which  is  silt  size.   Another 
possibility  is  that  shales  are  mostly  made  up  of  silt  size  particles ^ 
or  still  another  reason  might  be  that  the  silt  size  particle  is  much 
more  stable  and  therefore  resistant  against  degradation. 

It  is  interesting  to  note  the  relatively  small  values  of  coefficient 
of  variation  for  the  Atterberg  limits.   The  Atterberg  limits  can 
therefore  be  estimated  with  reasonable  confidence  if  a  few  test  results 
are  available.   Assuming  that  the  plastic  limit  and  liquid  limit  are 
normally  distributed  (which  is  reasonable  as  shown  in  Figure  7)  it  can 
be  said  that  68%  of  all  values  will  be  within  one  standard  deviation 
of  the  mean,  i.e.,  the  range  of  values  for  plastic  limit  is  19.9  to  24.1 
and  25.4  to  41.4  for  the  liquid  limit. 

It  is  clear  from  Figure  7  that  some  distinction  can  be  made  between 
low  and  high  durability  shales  based  on  the  histograms.   The  dispersion 
of  the  data  is  however  significant.   It  is  possible  that  further 
analysis  based  on  separate  geologic  groupings  will  distinguish  between 
strong  and  weak  or  durable  and  nondurable  shales. 

7. 2  Bivariate  Correlation  Analysis 

In  order  to  determine  whether  any  simple  correlations  can  be 
established  between  the  different  parameters,  bivariate  correlation 
analyses  were  performed  using  the  SCATTERGRAM  routine  of  the  SPSS 
package.   This  program  determine  the  statistics  indicated  in  Table  10 
as  well  as  a  two-dimensional  plot  of  data  points. 
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TABLE  10 


Statistics  Available  from  SCATTERGRAM 


Correlation  (r) 

The  Pearson  correlation  coefficient  r  is  used  to  measure  the 
strength  of  relationship  between  two  internal-level  variables. 
The  strength  of  the  relationship  indicates  the  goodness  of  fit  of 
a  linear  regression  line  to  the  data. 


r  = 


N         r   N   w  N 

y  x.  y.  -    y  x.      y  y. 
iii      ,**,  i    .**,  i 


i=l 


1=1 


i=l 


/N 


(r  N   9    f  N   ^|2  -If  N  „    f   N   -)2  .-])  1/ 


If  the  value  of  r  approaches  +  1.0  or  -  1.0,  we  can  assume  there  is 
a  strong  linear  relationship. 

2 
r  -  this  explains  the  proportion  of  variance  in  one  variable 

explained  by  the  other. 

Significance  tests  are  reported  for  each  coefficient  and  are 
derived  from  the  use  of  student's  t  with  N-2  degrees  of  freedom  for 
the  computed  quantity 


m 


in 


4.   Standard  error  of  estimate  -  this  statistic  is  the  standard  deviation 
of  the  residuals.   Residuals  are  the  errors  made  in  predicting  Y 
from  X  by  use  of  the  regression  equation. 


S  = 
e 


N 
1=1  1 


-   a 


N 

I   Y. 
i-1  ^ 


N      Nl/2 

b    y  X.  Y. 

i-1  X   X 


N-2 
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TABLE  10  (Continued) 


5.   Intercept  (a)  of  the  equation 


Y  =  a  +  b  X 


t   N     N  2^    f   S 


i=l   i=l 


1   X±   I  X  Y 
i=l   i=l 


2i      r  5     ^i2 


N  I      Xt 
i-1   X 


I  x. 

i-1 x 


6.   Slope  (b) 


b  = 


N  I  X  Y  -    I  X    J  Y 
l  j=l  x   ^    lj=l  x  j=l  x 

-   N    •>  r  N   -)2 
N  I   X  2  -   ^  X 

L  i-1  J  k-1  J 


7.  Plotted  values. 

8.  Excluded  values  (caused  by  changing  the  scale  of  the  plot) 

9.  Missing  values. 
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The  above  analysis  was  performed  for  all  combinations  of  the  varia- 
bles, and  plots  were  obtained.   The  values  of  r  (correlation)  and  the  number 
of  data  sets  with  reported  values  for  both  the  parameters  are  given  in  Table 
11.  The  plots  were  inspected,  and  it  was  decided  to  include  all  those  with 
r  >^  0.8,  even  those  which  involved  no  obvious  engineering  logic.   It  was 
further  decided  to  include  complete  data  on  all  the  analyses  which  gave  a 
'good'  correlation,  r  ^  0.7  was  used  as  a  cut-off  point.   These  results  are 
given  in  Table  12. 

The  compaction  data  are  all  generated  in  the  standard  AASHTO  test,  but 
the  maximum  size  was  in  some  cases  the  No.  k   sieve,  and  in  other  cases  the 
3A  in.  sieve.   The  data  were  examined  separately,  but  since  no  trend  was 
recognized,  they  were  lumped,  viz.,  Figures  15,  18,  20(b)  and  20(c). 

The  plotted  results  can  be  divided  into  three  important  parts,  i.e., 

a)  Good  straight  line  correlation  between  parameters  -  well  enough 
defined  over  a  wide  enough  range  so  that  results  of  the  one  can 
be  used  to  estimate  values  for  the  other.   This  can  be  seen  from 
Figures  9,  10,  12,  16,  17  and  18. 

b)  A  definite  trend  between  two  parameters  although  the  correlation 
is  low.   This  type  of  correlation  is  useful  in  determining  quali- 
tative relationships  which  may  assist  in  future  research.  An 
example  of  this  is  given  in  Figure  20. 

o)  After  plotting  the  data,  certain  boundaries  may  be  established 

within  which  all  the  data  presently  available  fall.   This  type  of 
diagram  may  be  used  in  the  future  to  determine  the  reliability  of 
laboratory  testing.   The  boundaries  may  change  when  more  data  be- 
come available  and  new  'groupings'  may  be  identified.   Examples 
of  this  type  of  plot  are  given  in  Figures  8  and  21. 


£ 

4> 


O 

m 


3 


4) 

■p 


to 
■p 

4! 
CO 

•) 


' 

a 

M  o 

mm 

<* 

g 

„8~2 

c 

w\«>  «\  r* 

o     o. 

•*   —  -*    —  *"" 

•4 

ipci  m  o  */*v 
©      *-•      n 

8 

009" 

«rs                  «s        jj 

ro        n        »—       r- 

g 

n  m  «9  c\  <o  ia  a  0 

j»  ftf  r>nwno 

O       0       p       p 

1         l         1 

(ft          W          O          «N          1 

1)         (O          - f          Ch         O 

«©^if»_3>r'*mn"»o<*io 

£ 

_o  ONO«*>tr\r-»r»(^»^-» 

?'?'?'?    ° 

«5       «       _j        ia       ^       0 
«  **i  t\i  cn  cj  0  «—  rtr-  » j  r-  m 

n 

_»  C4  -s  vn  ^»  i/\nonj  ^»  *\t 

p-        «j        p        O        d        p* 

j»              j»      e*     0      *- 

tr%       j*        n\        *-t       0        pj       «* 

s 

OO  JJ  w  J   •-*  -a   O  «-»  cj  n  0  r->  O 

OOOOOOO 

2 


§ 

8 


■-«        W        *—       -r         ©»        «-i        c*        «r\ 

j»        <0        f-         o        r—        r-       *o        r- 

60  *?    d    *?'    °    °    °    °    •? 

SO  V\  tf\  tf\  -J3  «J 

*C>  *2  «0   O  oj  O  -•   ;3  -*   2  "*  O  f\<D  r"i_»   fi  -* 
OOOOOOpOO 

Ar  ***  ^  ^  -*  u-\  <r>  so  ■-.  1— 
_**>-*  «~>  c\t  o  o  <o  (—  (O  NO 
«Hfir-'-»ih-^vD  r-  »-  o  1—  -*  r-Cs<o  00  00  cr\ 
WOOjiOhido  joj    ~-t  _r   O  fi  „  *-i  o  r-»  ^ 

oopooooooo 
•       It*       I 

S-*  VO  ^-J»  ^T  CD  li"\  NO  «D  _^  VO 

^b)h«o  ncod  ^vj  j  ioj  o^  ("  nvo  w  ca  n  pi 

*•         »-4  O  #H  Ci  W\  O  --t  a  .  - 

ooooooopopo 

.  O*  _  -•    Ok  #h    w~i    _j    rj     1 — 1    r*\    \n    rti    o— . 

•-t  •«     »^OOfN*^,C>^4«^m^_, 

opoQoopddood 
■   ■       1   1   *   l  111 

JO  **-*■)  CM         \0  CM         _s-         «n  ^1  u-\  m 

trgw       2       so        «-t       <q        i^r-eovJ^tvao 
«OO^tWO<H^*c>v^iC>.^^rUtrt^   r-i_»   u-»  VO  -1  NO  .X  SO   M 


r»       it*  o\      vo       i—       *-*       _^        ■»• 

•h         o  o        _»         0\.        _J-         Of—. 

CM  -»  «  o  to  rj  1 1-  f«<  i^r-  rih-  <m  no  eo  f— 

■0  '"'WJ  -«  flo  w^O  Oi^Q  «->eO  w-veo  ^  »»  < 


—  CMCM^-Ovj-J^oSn.'Xj^'SS  ?"*         T*i 

ft*OOo\iAtort»w»oioooooja)f4>OH4  hvo^^w^o^j- 

op'p'oddoo'dpddddd 


°   ° 


SfM       no        <d        _*        r*       ay        ^        j        ^ 
oooooddddd 

!i?-*r?_°       *°        **y        *"*        **•       <o        to        u-i        ri 

T?  ■?  <?  ■?  i  i  <?' 


u 

a) 

1 

as 

41 

.c 
•p 

§ 
§ 

■H 

4) 

o 
a 


<<frt<f',°        O        «j        o        d        tj 


«o        WN       *4        r-       ,rf         cn       c 

•-»         P—        rl         o         o         ^t         . 
^ri^r-Or-r^vOrg   (MO  f(M   C*cori_»    r-»\ 


'**>       o\      «j       cj       v\       ot       r~       r—       r-       rv 

«.  S^tT        ^        ""*        f~        c\        ^        cm        n        ^ 
i^CAH<OPjWrtCAHioo<on(DH'Ooij  jxjoirtr*^o»oH 

"c*"°"o"p"o   ddodddp'ddddddod 

eoop'doo^pop'oodddddp'dd 

oop*p*oodp'd<j*^d^o"pdp*pdddd 


•J^   fg    (M 


OcOOOOcOOnOC^OOtOOOr—  OJ-^ 


«?      <? 


•-^w^j        o        ^-t        irv        tr»        n\o         cw         o         n~i        »-_» 

Qopoooop'ddop 


C 

a> 

C 

*H 
Hh 
%-* 
<U 
O 
O 


«4_  m      SO       K       Cv      <N       w% 

»tOajco«co-j  ^^-w  ON.J  ^h.j  c\._ 

CD  ^"  ^O  CJ  m  'j-v  .  i 


Oj  J  CD  JCO^O  Cj»0  CJ*nCTv-»<a  fn 


NO  V\  p4 


hhco       e\i\o       (\iVOso        »noo_»        ^t  ^  o  •-« -»       o       H       5\      \o       w 
oooppdoddddddodododd 


OOi/NaSOvjrsiOCVif— 04 

O  H  »^v        so  f 


H>«  MCOCO  J. 


?     «?     ? 


O  f— 

\  i->  m  w-v 

mnnn 


CO 

§ 

H 

6) 

> 


*-S       HN    (VC 

•<     f  iir-l 


t>       o  1 

•OlNH. 

>  t-  f-  r-  ^sc 

W    ° 

>  CO         1 
*        Cu         * 

I  O  ( 


•    O    rt    —    -. 


£  a 


»A      M       p       «       U       B      R 

cy      c       a,       &.       0,       cl.       5 


1  -T   J*    C7\_»     , 


>       o       O       1 

Ml 

ft    2    S 


2  e 

6    6    3 


61 


0) 


o 

Oh 


m  co  m  co  o  co  vo 
r^   r-»   r-   cm    v£)   r-   r-J   O 


i-\   O   iH   .H   «H 


vro>N(NooconNroNconcoo 
rsi        ojogoooMXi^cooinvrin 


co  <r    co  m  -<r      cm      ts  cm  \o  rs  m  .h  cm  oo    ro  r-      vo  co 

^co^i^Ovo^ovOvor^Nino^coonvornvoco-jnvoninNcsiin^^^inrH 

OOOO-— (OOOOOOOOOO^TrHinroOcO^D^OCTiOOO^JOf-lOOOtH 
OOOOOOOOOOOO^HOOO^CvlC^rg^HOOOOOOOOOOOOO 


0) 


n  co  h  m    n  m  o\  vo  ro  vo  vo  co  >n    vj  <j-  co  co  r^  m 

ninn<tinrsHNfino>o>co>jcootNj<tcnrjc>inno>-Tr-jHrinHOHcoc7\n 
coccNc>oJlnn•JcOH^OHHO^lnroc^o^^r0^oofvJco^^coo^c^vIHlnco^oc^ 
^ovcinc^c^c^c^vOxONcoinNNrco^o^fnr^vDvtvD^nH'jJOin^ONco^o 


oooooooooooo 
I  I  I  I  I  I  I 


o  o 
I 


co  r>-  r^  r^  r-* 

I  I       t-t       |  |       r-i 


O    CM    CN 
I 


O-O.HOCNO-HOOO'-H 


in  n  cT>  m  cj\  rv  on  vo  r-*.cocococNCX\^ocrtin 
ONMnoONnmcocsruorON^vjNCNi 
r-.  co  o>  <r   a\cOrHvrOsrvooi\ONCnrNfsjM3H 


vt    r>.   vo   H   H 

o  o  -j-  <r  ^o  -h 


cN)sorov3-r^.ro<rco<T<r<rcNcNco<rin<r 


CN   co    O   O    n   (N 
H    sf    CN    CM    00    H 


CO  CN 

mCT\CncOO>!-3-CN--3-CO-v3" 

ocoNfcoc^cnt^r^vo^ 
<j-cooOvOco,-ho,n 


rHrHOiHrH^HiHinOoO 


o 

HI 

M 

u 

U 

n 

a 

**-* 

E 

C 

Ul 

01 

• 

to 

-a 

h 

o 

a 

nt 

Vj 

V 

> 

u 

a. 

m 

0) 

» — 

a 

sg 

-a   o 

C   TJ 

a. 

0)      • 

C     CO 

i-t     > 


mvo^aininoHnNo>invro> 
vDmr^-^-^nNinrNvOvocorsrs 
m^SvrrHcoHHN^iocoroin 


<r  o^  co  m  cm  vTcoo>roino><tcoNrno 

oinro-<rovoco^r-inr--nHrOi— t  i-h  cm  *-o  co  on 

^nvtiNrsrHnmvoOHinco^aMn^oJH 


I 


Nxjinmc^cNN^ino^coinocooHinco 

vOuO-HOO^C^COcOMCNuOrHO^CO<f-CNO»HvO 
r-{     r~\  rH  -H       |       CM     ,H       I 


NvJOOOCM!MsjNTtM<IONin 
COCO    HH    CM    (M  sf    H    H 


o>  <r  c~*  co  m  m  co         r^coinvOOic^vtHcotNco 

nnMno^M^tM^vOr- (.— 40v0r-nr--cor-.coo 

incoono^oninmcocooomconcoroMsrco 


rs  o^  \o  <t  Is- 
ro  r-*.  f~\  f-n  cm 
vT  un  o~>  vo  <r 


cm  rO  ^O  rO  vO  uO  CM 
iH  r-  in  On  cr*  r-.  o 
VO  <T   CM   O    r*»   CO   VO 


r-.  C"\   .—)  cr> 


-<r  o  co  cm 

CM    CM    t-(    CM 


r^  o>  O  CT\  cm  cr>  r- 


<r-jcocMO^ococT\r--(Mcoinin 

CNCN.HCOCOCNCNCNCN 


<■  <r  nnrs  vj 


vO 

i— i 

vo 

O 

t  3 

o\ 

O 

ON 

o 

CM 

H 

r*1  fO  ^O 

r^- 

r^ 

r^ 

r~  oo  oo  co 

r^cscocOi— tcorsi  in  in  (n^ 
^ai^noiOMntoontooi 
co^r^rgo^tHc\ioor-~r-HCNiinvo 


m 
ficomcomco<roominm 

Ol^t— (vOr-Hr^rH^OCT\00(N 

CNWOfflHOirJMO^O 
N\O(JiiOrN^Is00N(^C0 


oooooooooooooooooooooooooooooooooo 
I     I     I     I     I     I     I  I     I    I         I     I  II  I 


MfNincNNiniNinmNNin        m  o 

aQQQQQaQQQQQSQWl-ISSJ 


—i  ai  as  o       o  ps.-toa 

Oi3UUO'CXiCl«c/JCuU5QOQO 


rlHHHNNNNQnQQOM»101«l 

i-HhHh-jH- it— immi— iMNininmNminin 
oinuioiMoiuiMBiMinuiwoioioiwXzgsxSE 


J    ,-J    J    J    l-l 


i  £ 

a, 


^ 


PQ 


6? 


„  O 


o 
o 


<x> 


(IIQS)    I    310^0     X9pU|     6U()|D|S 


4 


•  I 


ro 

00 

o 

II 

k_ 

CVJ 

1—*, 

3 

CO 

00 

in 

a> 

(/) 

o 

> 

o 

OJ 

-a 

c 

Q 

co 

C*> 

C 

ro 
ro 

-a: 

a 

00 

LLi 

or 

Z> 

63 


• 

• 

• 

• 
• 

• 

• 
• 

•  • 

•  • 

•  •                   • 

• 
•                          •                   • 

• 

•                                      • 
• 

•           « 

•  • 

•  •   • 
• 
• 

• 

•               •       • 

• 
• 

• 

•               • 
• 

•                                                  • 

•                                                   • 
• 

•                  • 

• 

• 

• 

•  • 

•              # 

• 

•               • 

• 

• 
« 

t        •                                                  •                • 

• 

• 

• 
• 

• 

•  • 
• 

• 

H : : ' 

• 
i 

• 
• 

.-■ 

^r 

Q 

GO 

o 

CO 

■ — 

II 

^ 

Q 

CM 

> 



0) 

a 

cr 

(0 

o 

o 

ST) 

en 
> 

X 

V 

C 

Q 

IT) 

n 

>N 

CD 

!c 

o 

u 

ro 

3 

• 

rO 

Q 

CD 

UJ 

D 

3 

CO 

o 
o 


r^ 


ro 

ro 


(2IQS)  S  9P^D    xapui    6u|>|D|S 


6k 


•  • 


o 
o 


CO 
CM 

Q 

CO 


o 
h-      O 

> 

a: 


o 

o 

<M 


:-0 


O 

o 


1^- 

10 


fO 
fO 


(2IQS)  S  9P^0    xapui    6u!>jD|S 


00 

o 

i 

ii 


OJ 

Q 
CO 


00 

> 


C 

CO 

Q 

<n 

>. 

*; 

15 

R    ? 

o 

fO       3 

— > 

Q 

LU 

01 

(Z 

J£ 

D 

C 

CD 

CO  — 


I  • .  • 

:  • 


r-2 


CO 

in 

Q 
CO 


-a 

O) 

f-    o 
to    CO 

> 

cc 


ro 
rO 


o 

o 


X 

o> 

c 


o 

k. 

Q 


CO 


L^-O 


o 

O 


CD 


ro 
ro 


(2IQS)  S  9P*0    *»PU|    6u!>|D|S 


66 


•I 


--O 


to 

hi 
Q 
CO 


■X3 

n:   o 

&    CO 


> 
cr 

o 

o 

IT) 

x 

<y 


ro 


o 

O 


CO 


O 
O 


CO 


ro 
ro 


(QSQS)  *JQ   Aay  009   *9pu|    ^illiqojna    9>|D|S 


•  -o 


67 


1    .• 

• 

• 
•                    • 

• 
• 

•      a 

• 

• 

• 
• 

• 

• 
• 
•      • 

• 
• 

t 

• 
• 
i 

9 
• 

* 
• 

• 

• 
• 

H 

• 

• 

• 

• 

* 
• 

• 
• 

• 

•                • 

• 

• 

•                                                                              • 

• 

• 

• 

• 

• 
ft 

•                                                                      •               • 

•      • 

•                          • 

L± 

o 
o 


CO 

ID 

^^ 

Q 



CO 

CO 

" — 

o 

-a 

a> 

II 

JC 

k_ 

K 

a 
o 

" 

U> 

CO 

CO 

> 
■1) 

Q 

en 

(/) 

o 

o 

CO 

m 

> 

X 

<v 

CO 

■a 

ID 

a 

>> 

CO 

.t: 

J5 

PO 

o 

ro 

to 

Z3 

Q 

UJ 

<w 

tE 

D 
(0 

CO 

o 
o 


fO 


(S2QS)  P9>|dos  Aay  002   xepu|   AunqDjna    e>|D|S 


68 


CO 


GO 


en    _ 
o 


O 

o 


CO 


ro 
to 


(SW)  sson    DJ9d   *|S9i     sssupunos    pa^ipo^ 


...  o 


6q 


•  • 


o 


00 

— 

O 

Q 

ii 

(£> 

£ 

L. 

C\J 

* — 

U_ 

<n 

a. 

2 

- 

en 

>» 

> 

^~ 

to 

C 

— — 

0> 

n 

O 

£ 

<_ 

<D 

£ 

m 

E 

OJ 

3 

E 

LU 

— 

or 

X 

D 

Z> 

o 


o 
o 


-4 
o 


00 


(SW)  ssoi    0J9d   '4S91     sseupunos    pai^ipo^ 


TO 


E 


o 


(f) 

> 


CD 


D 


(Id)  X9PU|    AiPUSDId 


71 


•  * 

t     • 


£ 

3 


a; 

D 

+ 

Q 
.   Q 

<X)    ii 

II 


o 
m 


CD 
r0 


00 

CD 

6 


9     11 


<*■  t: 


Q 
Q 


W    2 

—    o 
2 


4-o 


C\J 


(QM)    Pd   «  Ajisuaa    43M    IDJnpN 


72 


•        • 

• 

•   •   • 

• 

• 
• 

•  • 

• 
• 

• 

•  • 
•    • 

• 

\  • 

m 

•  • 

•  • 

•  • 

• 

• 

c 

0) 

c 
o 
<J 

L. 
OJ 

+ 
5 

Q 
QJ  JL 

11 

•         • 

• 

•    • 

• 
• 

• 

• 
• 

J . 

• 

1- 

| 

(0 

CM 


<fr 

<7> 

O 

^— ^. 

^H 

ii 

Q 

L. 

«■ 

O 

"•-* 

*»^ 

... 

<♦- 

o 

o 

CL 

£ 

- 

(/■) 

>> 

> 

^~ 

c/> 

c 



0) 

n 

o 

Q 

>N 

h. 

Q 

00 

£ 

CM 

o 

3 

e 

u 
or 

o 

3 

o 

CD 


O 


CD 
CM 


CM 


00 
0> 


(IQM)   pd    '  ^4!SU9Q    |9M    wnwjXDfli 


73 


•  ••    • 


c\J 


--O 


CM 

fO 


CVJ 


C\J 


C\J 


(IU90)   0J9d    '  U80    papDdiuoo    s\/ 


CD 

o 

CO 

II 

IX. 

CD 

(  1 

— 

cr 

CD 

o 

O 

a> 

a. 

CO 

- 

> 

CE 

CD 

CVJ 

O 

tr 

CD 

en 

C 

O 

.*: 

D 

O 

C/5 

0*> 

<p 

UJ 

cr 

< 

Z> 

o 

V* 


■i-s- 


a_ro 
ro 


O 
O 


fO 

ro 


-f 
O 


CM     c 

CVJ        Q_ 


Q- 


K 

c 


O 


--o 


CO 

o 

I 

H 


Q 
If) 
Q 
</) 

</) 

> 


t;     £ 


o 

CM 

UJ 


(QSaS)  *JQ  Aay  00S  xepui   A*mqDjna    9)|D|S 


75 


- a.» 


00 

o 

I 

M 


m 


I   - 


m 
o 

in 

> 


o 

CVJ 

Id 
K 

3 


4- 

o 

2 


iO 


to 
to 


(SSQS)  P«|do$  A*y  00fi  x»pu|   Aimqojna    ©hois 


76 


1  =■ 

« 


<7> 


m 


«> 

w 
■D 

*- 

o 

E 
£ 

Q. 

o 


CD 


Q 

OJ 
Q 
if) 

U~) 
> 


o 

CNJ 
UJ 

cr 
•=> 
o 


o 
o 


ro 


(Q2QS)  **Q  *»a  002  *apu|    ^IHqojna    a>|D|S 


77 


•  • 


q_c\j 


ro 

* 

O 

i 

u 

X 



00     (f) 

Q 

(\J 

<_> 

n 

0) 

<n 

Q- 

in 

<U 

> 

o> 

o 

I 

c: 

CD 

k. 

£ 

(O 

TJ 

O 

<»_* 

c 

o 

*  -J 

CJ 

LU 

q: 

z> 

o 

-■o 


o 
o 


ro 
ro 


(Q2QS)  *JQ  A»d  002    X9PUI   **!l!qDJna    9>1DIS 


78 

From  the  summary  in  Table  11  it  is  clear  that  very  few  good  corre- 
lations exist  on  a  one  to  one  basis  for  all  the  shale  data.   It  should, 
however,  be  remembered  that  the  Indiana  shales  cover  a  wide  range  of 
characteristics.   With  this  in  mind  it  was  decided  to  analyze  some  of  the 
data  after  dividing  the  data  sets  according  to  the  geology. 

From  the  description  of  the  geology  the  following  division  was  used: 

a.  Ordovician  -  all  formations 

b.  Devonian  -  all  formations  plus  riew  Albany  of  "ississippian 

c.  "ississippian  -  Borden  ^roup 

Locust  Point 
New  Providence 

d.  'Tississippian  -  rest  of  formations 

e.  Pennsylvanian  -  all  formations. 

Hot  all  the  parameters  were  analyzed  but  only  those  containing  some 
notion  of  durability  and  strength.   Table  13  indicates  which  analyses  were 
performed. 

Home  results  included  very  few  data  points,  and  the  results  were  con- 
sidered to  be  significant  only  if  more  than  9  data  points  were  included 

2 

and  the  value  of  r   >  O.65.   This  division  is  very  arbitrary,  but  it  is 

felt  that  values  obtained  by  this  method  can  surely  be  considered  as 
'significant'.   These  values  are  given  in  Table  1*+.   The  results  from  this 
Table  indicate  that  there  exists  reasonably  good  correlation  between  the 
different  slake  durability  measures.  When  these  results  are  compared  to 
that  in  Table  11,  it  is  clear  that  there  are  some  differences  between  the 
different  geological  formations.   The  basic  trends  are,  however,  still 
recognizable . 
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TABLE  13 

Combinations  of  Parameters  Analyzed  After 
Dividing  the  Data  According  to 
Geological  Differences 
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TABLE  14 


Significant  Correlations  of  Analyses  After 

Dividing  the  Data  According  to 

Geological  Differences 


Dependent 

Independent 

2 

r 

No.    of 

* 
a 

b 

var. (y) 

var. (x) 

Points 

a.    Ordovician 

SD5S 

SDI2 

0.714 

9 

99.0208 

-1.229 

b .    Devonian 

plus   New  Albany  of  Mississippian 

SDI2 

SDH 

0.979 

10 

0.1778 

0.1043 

SD2D 

SDH 

0.936 

9 

16.244 

-0.1602 

SD5D 

SDH 

0.9426 

10 

17.61 

-0.1754 

SD2S 

SDH 

0.946 

9 

13.32 

-0.1328 

SD5S 

SDH 

0.902 

10 

13.465 

-0.1337 

SD2D 

SDI2 

0.9865 

9 

156.943 

-1.5693 

SD5D 

SDI2 

0.8733 

10 

159.698 

-1.602 

SD2S 

SDI2 

0.9657 

9 

126.457 

-1.28 

SD5S 

SDI2 

0.8141 

10 

120.374 

-1.206 

SD5D 

SD2D 

0.788 

9 

3.6552 

0.9621 

SD2S 

SD2D 

0.9552 

9 

20.323 

0.806 

SD5S 

SD2D 

0.7166 

9 

28.007 

0.716 

SD2S 

SD5D 

0.875 

9 

28.707 

0.712 

SD5S 

SD5D 

0.9254 

12 

28.215 

0.722 

SD5S 

SD2S 

0.821 

9 

5.819 

0.929 

c.   Mississippian  -   Borden 

Group,    Locust   Point, 

New  Prov. 

SD5S 

SD5D 

0.6871 

24 

66.688 

0.292 

CBR1 

SD5D 

0.681 

16 

73.853 

0.733 

CBR1 

SD5S 

0.642 

16 

36.22 

1.9287 

d.    Rest   of 

Mississippian 

SDI2 

SDH 

0.6566 

61 

-13.484 

0.787 

SD5D 

SDI2 

0.6932 

52 

100.623 

-0.908 

SD2S 

SDI2 

0.6603 

36 

94.97 

-0.86 

CBR2 

SDI2 

0.697 

17 

104.26 

-9.552 

DD1 

SD2D 

0.656 

17 

-391.67 

3.788 

CBR2 

SD5D 

0.678 

18 

11.9599 

7.196 

CBR2 

SD2S 

0.729 

16 

4.84 

7.99 

CBR2 

SD5S 

0.701 

18 

-0.656 

7.126 

e.    Pennsylvanian 

SD5S 

SDI2 

0.631 

27 

92.62 

-0.957 

SD2S 

SD2D 

0.781 

13 

16.99 

1.092 

y  =  a  +  bx 
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The  good  correlations  between  the  different  slake  durability  measures 
are  to  be  expected  and  there  will  be  an  elaboration  on  this  aspect  in 
the  next  section. 

Inspection  of  the  different  plots  indicated  that  'boundaries'  can 
also  be  established  for  certain  parameters  after  dividing  the  data 
according  to  geology.   An  example  of  this  is  given  in  Figure  21.   This 
series  of  plots  of  Slake  Durability  Index  500  rev  dry  vs.  Slaking  Index 
Cycle  1  show  distinct  groupings  of  the  data  for  the  different  geologic 
formations.   This  presentation  not  only  gives  insight  in  the  typical 
properties  of  the  different  formations,  but  may  also  be  of  help  in 
determining  the  reliability  of  laboratory  results.   For  example,  if  a 
result  of  the  formations  in  Figure  21(h)  gives:  SDI1  =  90,  SD5D  =  90, 
then  there  is  a  high  possibility  that  this  answer  is  wrong.   The 
reliability  of  this  will  improve  as  more  data  become  available.   Also  the 
data  in  Figure  21(g)  may  be  broken  up  further,  which  may  explain  the 
single  point  in  the  upper  right  hand  corner.   This  point  may,  however, 
also  be  the  result  of  a  'faulty'  laboratory  test.   There  exist  therefore 
many  possible  applications  for  this  type  of  analysis,  and  it  should  be 
pursued  in  the  future. 

7. 3  Multiple  Regression  Analysis 

The  REGRESSION  routine  of  the  SPSS  package  was  used  to  carry  out  the 
multiple  regression  analysis. 

With  the  large  number  of  parameters  available  in  the  data  on  shales, 
various  regression  models  might  be  possible  to  predict  various  parameters 
using  the  others.   For  this  study,  attempts  were  made  to  predict  California 
Bearing  Ratio  (CBR)  from  some  of  the  other  parameters.   The  reason  for  this 
is  that  CBR  is  considered  to  be  the  most  time  consuming  test. 
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(c)   Mississippian-Clore,  Palestine  Formations 


FIGURE  21.   Slake  Durability  Index  500  Rev.  Dry  vs.  Slaking  Index  Cycle  1 
for  Different  Geological  Formations. 
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(f)   Mississippian-Haney,  Big  Clifty,  Elwren,  Sample,  Bethel  For»ations 

FIGURE  21.  (continued) 
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FIGURE  21.  (continued) 
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It  was  first  of  all  necessary  to  pick  out  those  data  nets  with  data 
noints  for  all  the  different  parameters  to  be  used  in  the  analyses.   Event- 
ually 2h   sets  of  data  were  used. 

In  order  to  determine  the  first  models,  a  table  similar  to  Table  11 
was  inspected  and  five  parameters  including  2  durability  measures,  Atter- 
berg  limits,  grain  size  analysis  and  one  density  measure,  with  the  highest 
one-on-one  correlation  with  CRRI  and  CRR2 ,  were  chosen.   The  models  were 
then  set  up  to  include  first  order,  second  order  and  combination  terms. 
A  stepwise  regression  was  performed  in  which  the  variable  with  the  highest 
correlation  is  introduced  into  the  model  first,  followed  by  the  second 
highest,  etc.,  until  the  standards  set  by  the  program  are  reached.   See 
Table  15  for  the  parameters  used  in  the  first  models  as  well  as  the  results. 
It  can  be  seen  that  one  parameter  (SD2C)  (DDl)  was  removed  during  the  step- 
wise regression  for  CBR1.   Also  it  can  be  seen  that  in  order  to  predict 
CiiRl  and  CBR2  ,  values  of  all  five  parameters  must  be  known.   The  testing 
to  determine  these  might  take  as  long  as  the  testing  for  CRR1  and  CBR2  , 
which  makes  the  above  analysis  really  an  academic  exercise.   Even  though 
these  models  gave  satisfactory  correlations,  it  should  be  pointed  out  that 
only  2U  data  sets  were  used  and  they  need  further  validating  as  additional 
data  sets  become  available. 

Another  model  was  set  up,  for  which  more  data  sets  had  the  values  of 
all  the  parameters  available.   The  original  2U  data  sets  were  used  to  set 
up  the  model  and  it  was  then  tested  against  the  data  of  Deo  (1972),  i.e.  , 
laboratory  numbers  1-15.   This  model  to  predict  CBR1  is  given  in  Table 
16.   Figures  22  and  23  indicate  how  well  the  model  fits  the  data  of  the 
2U  data  sets  and  also  the  data  by  Deo  (1972).   The  latter  is  over  a 
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TABLE  15 


Regression  Models  to  Predict  CBR1  and  CBR2  Based  on 
High  One-tc-One  Correlation 


Parameters  used  for: 

CBR1:  SD2S,  FN,  PL,  PCOL,  DDl 
CBR2:  SD5S,  FN,  PI,  PCOL,  PM 

Results  for  CBR1  (stepwise): 


Variable 

Multiple 

2 

r 

2 

r 

Overall 

Signifi- 

r 

change 

F 

cance 

(SD2S)(PC0L) 

0.3513 

0.1234 

0.1234 

3.097 

0.092 

(PCOL)2 

0.3780 

0.1429 

0.0195 

1.7508 

0.198 

(PCOL) 

0.4638 

0.2151 

0.0722 

1.8268 

0.175 

(SD2S)(DD1) 

(SD2S)2 

(DDl)2 

0.5267 

0.2774 

0.0624 

1.8238 

0.166 

0.582 

0.3388 

0.0613 

1.8442 

0.155 

0.6306 

0.3977 

0.0589 

1.8708 

0.145 

(SD2S)(PL) 

0.6759 

0.4568 

0.0591 

1.922 

0.132 

(PL) (PCOL) 

0.7137 

0.5094 

0.0526 

1.9465 

0.127 

(SD2S) 

0.7707 

0.594 

0.0846 

2.2756 

0.081 

(SD2S)(FN) 

0.7918 

0.627 

0.033 

2.1848 

0.094 

(FN) 

0.7983 

0.6373 

0.0104 

1.9169 

0.139 

(FN) (PL) 

0.8789 

0.7724 

0.1351 

3.1116 

0.035 

Remove 

(SD2S)(DD1) 

0.8788 

0.7723 

-0.0001 

3.7009 

0.017 

Add 

(FN)2 

0.9088 

0.8259 

0.0536 

4.3485 

0.01 

(FN) (PCOL) 

0.922 

0.8505 

0.0254 

4.3743 

0.012 

Final  model  with  r  =  0.922: 


* 

* 

Variable 

B 

Variable 

B 

(SD2S)(PCOL) 

0.001553 

(SD2S) 

0.2382 

(PCOL)2 

0.02454 

(SD2S)(FN) 

-0.004075 

(PCOL) 

-2.5188 

(FN) 

0.9224 

(SD2S)2 
(DDl)2 

-0.0028 

(FN) (PL) 
(FN)2 

-0.03936 

0.0009553 

0.002976 

(SD2S)(PL) 

0.001369 

(FN) (PCOL) 

-0.004183 

(PL) (PCOL) 

0.06984 

Constant 

1.3565 
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TABLE  15  (Continued) 


Results  for  CBR2 

(stepwise) : 

Variable 

Multiple 

2 

r 

2 
r 

Overall 

Signifi- 

r 

change 

F 

cance 

(SD5S)2 

0.6441 

0.4149 

0.4149 

15.6 

0.001 

(FN)(PM) 

0.6806 

0.4632 

0.0483 

9.06 

0.001 

(SD5S)(PC0L) 

0.6916 

0.4783 

0.0151 

6.112 

0.004 

(SD5S)(PM) 

0.7189 

0.5168 

0.0385 

5.08 

0.006 

(PCOL) 

0.7476 

0.5589 

0.0421 

4.561 

0.007 

(PCOL)* 

0.7676 

0.5892 

0.0303 

4.0636 

0.01 

(FN) (PI) 

0.7923 

0.6278 

0.0386 

3.8549 

0.012 

(PI) 

0.7965 

0.6345 

0.0067 

3.2543 

0.023 

(FN) (PCOL) 

0.8097 

0.6556 

0.0211 

2.9609 

0.034 

(PI) (PCOL) 

0.8158 

0.6656 

0.01 

2.5874 

0.056 

(FN)2 

0.8184 

0.6698 

0.004 

2.2131 

0.094 

(FN) 

0.8363 

0.6995 

0.0296 

2.1334 

0.110 

Final  model  with  r  =  0.8363: 


Variable 


(SD5S) 

(FN)(PM) 

(SD5S)(PC0L) 

(SD5S)(PM) 

(PCOL) 

(PCOL)* 

(FN) (PI) 

(PI) 

(FN) (PCOL) 

(PI) (PCOL) 

(FN)2 

(FN) 

Constant 


B* 


0.0006265 
-0.005126 

0.008285 
-0.01236 
-0.6047 

0.005335 

0.024 
-1.4375 
-0.01171 

0.03876 

0.001742 
-0.2141 
23.6259 


*  CBR2  =  B  (SD55)' 


B  (FN)(PM) 


B  (SD55)(PCOL)  + 
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TABLE  16 


Model  to  Predict  CBR1,  Set  Up  for  24  Data  Sets  and 
Tested  for  15  Different  Sets 


1.  Parameters  used: 

SD5D,  FN,  PI,  PSI,  PM 

2.  Results  for  CBR1  (stepwise) 


Variable 


Multiple 

r 


r 
change 


Overall 

F 


Signifi- 
cance 


(SD5D)(FN) 

0.2816 

0.0793 

0.0793 

1.8952 

0.182 

(PM) 

0.4109 

0.1688 

0.0895 

2.1326 

0.143 

(FN) (PI) 

0.4593 

0.211 

0.0421 

1.7823 

0.183 

(FN)2 

0.5207 

0.2711 

0.0602 

1.7666 

0.177 

(FN) 

0.6808 

0.4636 

0.1925 

3.1107 

0.034 

(PM)?, 
(PI) 

0.7256 

0.5264 

0.0629 

3.1495 

0.029 

0.7343 

0.5392 

0.0128 

2.6746 

0.049 

(PD(PM) 

0.7795 

0.6076 

0.0684 

2.9034 

0.036 

Remove 

(FN) (PI) 

0.7795 

0.6076 

-0.00001 

3.5393 

0.017 

Add 
(SD5D) 

0.7959 

0.6335 

0.0259 

3.2411 

0.024 

(SD5D)(PM) 

0.8129 

0.6607 

0.0272 

3.0296 

0.031 

(PI) 

0.8317 

0.6918 

0.031 

2.9178 

0.037 

(PI) (PSI) 

0.8647 

0.7477 

0.0559 

3.2324 

0.028 

Remove 

(PM) 

0.8647 

0.7476 

-0.00004 

3.8512 

0.013 

(psir 

0.8704 

0.7575 

0.0099 

3.408 

0.023 

(SD5D)(PSI) 

0.8734 

0.7629 

0.0054 

2.949 

0.042 

Remove 

(PI)2 

0.8734 

0.7628 

-0.0001 

3.5078 

0.020 

Add 

(FN) (PM) 

0.8807 

0.7757 

0.0129 

3.17 

0.033 

(FN) (PSI) 

0.8948 

0.8007 

0.0250 

3.0906 

0.04 

Final  model  with  r  =  0.8948: 


Variable 


Variable 


(SD5D)(PI) 

-0.005012 

(PI) 

-1.4378 

(FN)2 

-0.00431 

(PI) (PSI) 

0.01197 

(FN) 
(PM)Z 

0.50302 

(PSI)2 

-0.004016 

-0.004572 

(SD5D)(PSI) 

0.0036007 

(PD(PM) 

0.07577 

(FN) (PM) 

-0.01907 

(SD5D)2 

0.0009055 

(FN) (PSI) 

0.0031716 

(SD5D)(PM) 

-0.01654 

Constant 

11.927105 
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TABLE  16  (Continued) 


3.   Results  for  CBR2  (stepwise) 


Variable 

Multiple 

2 
r 

2 

r 

Overall 

Signifi- 

r 

change 

F 

cance 

(SDI2)(PH) 

0.7016 

0.4922 

0.4922 

21.324 

0.000 

(SDI2)2 

0.7503 

0.563 

0.0708 

13.5247 

0.000 

(PH) 

0.8237 

0.6784 

0.1155 

14.0647 

0.000 

(XY)2 

0.826 

0.6822 

0.0038 

10.1976 

0.000 

(SDI2)(XY) 

0.8417 

0.7084 

0.02619 

8.7462 

0.000 

(DD)2 

0.8513 

0.7248 

0.0164 

7.4612 

0.000 

(SDI2)(AMC) 

0.8594 

0.7386 

0.0138 

6.457 

0.001 

(DD) (AMC) 

0.8706 

0.7579 

0.0194 

5.8702 

0.002 

(XY) 

0.8918 

0.7954 

0.03745 

6.0459 

0.002 

(SDI2) 

0.9091 

0.8265 

0.0312 

6.1947 

0.002 

(PH) (AMC) 

0.9138 

0.8351 

0.0086 

5.5242 

0.003 

(PH)(XY) 

0.9175 

0.8418 

0.0067 

4.8759 

0.007 

Remove 

(DD)2 

0.9175 

0.8417 

-0.00001 

5.8022 

0.003 

Add 

(AMCr 

0.9186 

0.8438 

0.002 

4.9507 

0.006 

Remove 

(DD) (AMC) 

0.9186 

0.8438 

-0.0000 

5.8917 

0.002 

Add  9 
(PH)Z 

0.92 

0.8464 

0.0026 

5.0507 

0.006 

Final  model  with  r  =  0.92: 


Variable 


(SDI2)(PH) 

0.012821 

(SDI2)2 

0.0004225 

(PH), 

6.80998 

(XY)2 

-0.17319 

(SDI2)(XY) 

0.00738 

(SDI2)(AMC) 

0.02098 

(XY) 

4.2569 

(SDI2) 

-0.37786 

(PH) (AMC) 

-0.35246 

(PH)(XY) 

-0.30195 

(AMC)2 
(PH)2 

0.04786 
-0.22927 

Constant 

-21.7448 
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larger  range,  and  the  fit  (r  from  regression  analyses)  is  better  than 
that  of  the  model  although  the  regression  line  does  not  correspond  well 
with  the  line  x  =  y.   This  model  still  includes  difficult-to-determine 
parameters,  and  might  therefore  be  impractical.   It  should,  however,  be 
tested  as  more  data  become  available. 

A  final  model  based  on  simple-to-perform  tests  was  tried.   The  para- 
meters and  the  results  are  given  in  Table  IT.   Tt  can  be  seen  that  both 
models  are  very  good  for  the  2k   sets  of  data.   This  model  must  still  be 
validated,  and  it  is  proposed  that  this  be  done  as  more  data  become 
available. 

With  the  good  correlations  obtained  between  the  different  slake 
durability  parameters  on  a  one-to-one  basis  it  was  decided  to  analyze 
these  correlations  again  using  second  order  equations  and  the  REGRESSION 
program.   For  example,  the  relationship  between  SDI1  and  SDI2  would  be 
determined  by  the  equation: 

SDI2  =  C  (SDH)2  +  C  (SDI2)(SDI1)  +  C  (SDIl)  +  C^ 

The  results  of  the  different  analyses  are  given  in  Table  18.   It  is 
clear  that  the  correlations  are  much  higher  than  those  given  in  Table 
11  which  were  based  on  bivariate  analysis.   It  is  therefore  recommended 
that  these  second  order  equations  be  used  for  predicting  the  different 
slake  durability  parameters  when  one  is  known.   It  should  be  remembered 
that  these  models  are  based  on  all  163  sets  of  data  and  were  not 
validated.   New  data  that  become  available  should  be  used  for  validating 
these  models  so  that  their  reliability  can  be  increased. 


TABLE  17 

Regression  Models  to  Predict  CBR1  and  CBR2  from 
Simple-To-Perform  Tests 


1.  Parameters  used  (for  CBR1  and  CBR2) 

SDI2,  DD,  AMC,  PH,  XY 

2.  Results  for  CBR1  (stepwise): 


Variable 

Multiple 

2 
r 

2 
r 

Overall 

Signifi- 

r 

change 

F 

cance 

(SDI2)(PH) 

0.3909 

0.1528 

0.1528 

3.9672 

0.059 

(SDI2)(XY) 

0.6444 

0.4153 

0.2625 

7.4563 

0.004 

(DD)2 

0.6578 

0.4328 

0.0175 

5.086 

0.009 

(DD) (AMC) 

0.6951 

0.4832 

0.0504 

4.4409 

0.011 

(PH) 

0.7052 

0.4973 

0.0141 

3.561 

0.021 

(SDI2)(DD) 

0.7179 

0.5154 

0.0181 

3.0128 

0.034 

(PH)2 

0.766 

0.5868 

0.0714 

3.2458 

0.024 

(AMC) 

0.7757 

0.6016 

0.0149 

2.8318 

0.039 

(XY)(AMC) 

0.7787 

0.6063 

0.0047 

2.3959 

0.069 

(XY)2 

0.7944 

0.6311 

0.0248 

2.224 

0.089 

(DD)(XY) 

0.8099 

0.6559 

0.0248 

2.0793 

0.112 

(PH)(AMC) 

0.8167 

0.667 

0.0111 

1.8361 

0.162 

(SDI2)(AMC) 

0.8366 

0.7 

0.0329 

1.7943 

0.18 

(PH)(XY) 

0.8501 

0.7227 

0.0228 

1.6755 

0.22 

(AMC)2 

0.8555 

0.7319 

0.0092 

1.4561 

0.302 

(SDI2)2 

0.8657 

0.749 

0.0175 

1.308 

0.375 

Final  model 

with  r  =  0.8657: 

Variable 

B 

Variable 

B 

(SDI2)(PH) 

0. 

05855 

(XY)(AMC) 

-0 

.26375 

(SDI2)(XY) 

0. 

017099 

(XY)2 

0, 

.06404 

(DD)2 

0. 

0028604 

(DD)(XY) 

-0, 

.023976 

(DD) (AMC) 

-0. 

026762 

(PH) (AMC) 

-0, 

,25951 

(PH) 

-14. 

8706 

(SDI2)(AMC) 

-0, 

.01517 

(SDI2)(DD) 

-0. 

0030687 

(PH)(XY) 

0. 

,46754 

(PH)2 

0. 

84256 

(AMC)2 

-0, 

,142685 

(AMC) 

10. 

9822 

(SDI2)2 
Constant 

0. 
-6. 

0004864 
4060052 

9k 


3.   Results  for  CBR2  (stepwise) 


Variable 

Multiple 

2 
r 

2 

r 

Overall 

Signifi- 

r 

change 

R 

cance 

(SDI2)(PH) 
(SDI2) 

0.7016 

0.1+92 

0.1+92 

21.321+ 

0.000 

0.7503 

0.563 

0.0708 

13.525 

0.000 

(PH)_ 
(XY)2 

0.8237 

0.678 

0.1155 

11+.065 

0.000 

0.826 

0.682 

0.0038 

10.198 

0.000 

(SDI2)(XY) 

0.8U17 

0.708 

0.0262 

8.7^6 

0.000 

(DD)2 

0.8513 

0.721+8 

0.016I+ 

7.1+612 

0.000 

(SDI2)(AMC) 

0.859*+ 

0.7386 

0.0138 

6.1+57 

0.001 

(DD)(AMC) 

0.8706 

0.7579 

0.0191+ 

5.8702 

0.002 

XY 

0.8918 

0.795*+ 

0.0375 

6.0U59 

0.002 

SDI2 

0.9091 

O.8265 

0.0312 

6.191+7 

0.002 

(PH)(AMC) 

0.9138 

0.8351 

0.0086 

5.521+2 

0.003 

(PH)(XY) 

0.9175 

0. 81+18 

0.0067 

1+.8759 

0.007 

Remove 

(DDT 

0.9175 

0.81+17 

-0.00001 

5.8022 

0.003 

Add 
(AMC) 

0.9186 

0.81+38 

0.00203 

1+.9507 

0.006 

Remove 

(DD)(AMC) 

0.9186 

0. 81+38 

0.0 

5.8917 

0.002 

Add 
(PH)2 

0.92 

0.81+61+ 

0.0026 

5.0507 

0.006 

(XY)(AMC) 

0.9205 

0.81+71+ 

0.001 

1+.2712 

0.011+ 

(DD)(XY) 

0.9208 

0. 81+78 

0.0001+ 

3.5802 

0.030 

( DD ) ( AMC ) 

0.9222 

O.850U 

0.0027 

3.0328 

0.059 

(DD) 

0.9231 

0.852 

0.0016 

2.519 

0.110 

Final  model  with   r 

=  0.9231 

Variable 

B 

Variable 

B 

(SDI2)(PH) 

O.Oi+195 

(PH)(AMC) 

-0.57757 

(SDI2) 

0.005662 

(PH)(XY) 

-O.I+I+896 

(PH) 

/  it-**  \  2 

8.1201+6 

(AMC)2 
(PH)2 

0.086209 

{Ti)d 

-0.12581+ 

-0.26I+60I+ 

(SDI2)(XY) 

0.021707 

(XY)(AMC) 

-0.21+1+53 

(SDI2)(AMC) 

0.021051 

(ddMxy) 

-0.01+5907 

(XY) 

12.2991* 

(dd)(amc) 

0.018091* 

(SDI2) 

-0.691+905 

(DD) 
Constant 

0.13003 
-1+9. 9787 
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TABLE  18 


Results  of  Regression  on  Second  Order  Equations  for 
Durability  Parameters 


Dependent 

Independent 

Cl 

C2 

C3 

C4 

Multiple 

variable 

(y)  variable (x) 

r 

SDI2* 

SDH 

-0.03101 

0.03247 

0.36789 

19.1696 

0.9018 

SD2D 

SDH 

-0.01763 

0.01784 

-2.6874 

90.5737 

0.8984 

SD5D 

SDH 

0.02093 

0.01759 

-2.8427 

83.325 

0.8998 

SD2S 

SDH 

0.02414 

0.02177 

-3.0096 

78.699 

0.8981 

SD5S 

SDH 

0.02507 

0.0282 

-2.993 

69.224 

0.8399 

SD2D 

SDI2 

0.00614 

0.01205 

-1.6672 

97.144 

0.9798 

SD5D 

SDI2 

0.00956 

0.0142 

-1.97098 

92.5352 

0.9759 

SD2S 

SDI2 

0.01203 

0.01374 

-2.166 

93.406 

0.9808 

SD5S 

SDI2 

0.01478 

0.01638 

-2.3413 

85.7386 

0.9495 

SD5D 

SD2D 

-0.00888 

0.01223 

0.6633 

3.2595 

0.9899 

SD2S 

SD2D 

-0.00978 

0.01266 

0.7319 

2.063 

0.9808 

SD5S 

SD2D 

-0.00598 

0.01181 

0.4334 

3.901 

0.9732 

SD2S 

SD5D 

-0.00696 

0.0144 

0.1184 

20.87 

0.9045 

SD5S 

SD5D 

-0.00753 

0.01309 

0.4312 

7.8161 

0.9626 

SD5S 

SD2S 

-0.00914 

0.01236 

0.6928 

2.7344 

0.9795 

C,  x  +  C  xy  +  C.  x  +  C, 


For  example:  SDI2  =  (^(SDIl)   +  C2  (SDI2)  (SDH)  +  C3(SDI1)  +  C, 
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8.   Conclusions 

1.  The  storage  and  retrieval  system  described  in  this  report  is 
only  a  first   attempt     and  should  be  improved  in  the 
future  by  the  inclusion  of  additional  data  as  they  become 
available. 

2.  The  shales  of  Indiana  cover   a  wide  range  of  characteristics 
and  it  is  necessary  to  have  complete  data  sets  on  as  many 
different  shales  as  possible  in  order  to  obtain  good 
statistical  correlations. 

3.  The  values  of  all  the  slake  durability  indices  can  be  estimated 

if  one  is  known  by  using  the  second  order  equations  of 

2 
Table  18,  which  are  of  the  form,      SDI2  =  C,  (SDH)   + 

C0(SDI2)(SDI1)  +  C-(SDIl)  +  C. . 

/  J  4 

4.  Various  relations  exist  between  the  data  on  a  one-to-one 
basis  which  give  reasonable  correlations.   Details  of  the 
most  significant  relations  are  given  in  Table  12. 

5.  Various  regression  models  can  be  obtained  to  predict  some 
parameters  by  using  others.   Good  models  were  obtained  for 
estimating  CBR  by  using  5  parameters.   It  should  be  remembered 
that  the  data  used  for  the  models  come  from  all  the  different 
geological  series.   Better  models  may  result  in  the  future 
when  more  complete  data  sets  are  available  for  separate 
geological  formations. 

6.  This  partial  analysis  indicates  that  the  data  can  be  divided 
into  groups  based  on  their  geological  origin.   This  area 
should  be  expolored  more  fully  in  the  future. 
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9.   Recommendat  ions 

1.  A  general  sorting  program  should  be  written  so  that  the  sorting 
of  data  can  be  done  in  a  more  detailed  manner.   It  is  desirable 
to  have  a  program  to  sort  and  print  only  certain  data  sets  or 
parts  of  data  sets.  Sorting  of  this  kind  can  be  accomplished  at 
this  stage  by  using  some  routines  of  the  SPSS  package. 

2.  The  statistical  analyses  presented  in  this  report  should  be  seen 
as  first  approximations.   More  detailed  analyses  should  be  done 
at  a  later  date  when  more  data  are  available.   It  is  important 

to  create  models  with  some  of  the  data  available  and  test  it  with 
the  rest  otherwise  very  little  confidence  can  be  placed  on  the 
models  obtained. 

3.  The  notion  of  'grouping'  data  according  to  geological  origin  should 
be  further  explored.   It  is  necessary  to  obtain  more  complete  data 
sets  from  some  of  the  geological  formations  before  this  can  be 
done  with  high  statistical  confidence.   A  complete  project  in 

this  area  may  prove  to  be  one  of  the  most  useful  for  the 
practitioner.   Data  from  Illinois,  Kentucky  and  Ohio  may  help 
to  prove  the  findings  of  any  local  study. 

4.  The  maintenance  of  a  shale  data  bank  should  be  assumed  at  an 
early  data,  by  the  Indiana  State  Highway  Commission.   Periodic 
statistical  examination  of  data  from  this  bank  by  both  ISHC 
and  Purdue  University  is  envisioned. 
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APPENDIX  A 


DESCRIPTION  OF  TESINO  METHODS  USED  BY  ISHC 
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A-l.   Method  of  Test  for  Determining  Loss-On-Ignition* 
Test  Method  No.  Ind.  501-76 

1.  Scope 

1.1  This  method  covers  the  determination  of  organic  matter  (boss- 
On  -Ignition)  in  soils.   The  procedure  is  based  on  oxidizing 
all  of  the  organic  material  present  and  stabilizing  the  or- 
ganic residue  at  1200°F. 

2 .  Apparatus 

2 . 1  Electrical  Muffle  Furnace 

The  electric  muffle  furnace  should  have  good  air  circulation 
and  be  capable  of  having  its  temperature  regulated  at  1200°F. 
±50°F. 

2 . 2  Porcelain  Crucibles 

The  crucibles  shall  be  1-3/8  inch  in  depth  and  1-5/8  inch  in 
diameter.   Similar  dishes  of  platinum  or  silica  may  be  used. 

2 . 3  Balance 

A  balance  conforming  to  the  requirements  of  AASHTO  M231 ,  Class  C, 

2.1+  Oven 

A  thermostatically  controlled  drying  oven  capable  of  maintain- 
ing temperatures  of  110°C.  ±   5°C.  (230°F.  ±  9°F. )  for  drying 
soil  samples. 

2 . 5   Spatula 

A  spatula  or  pen   knife  having  a  blade  about  3  inches  (76  mm) 
in  length  and  about  3 A  inch  (19  mm)  in  width. 


*This  test  method  is  only  tentative  and  is  under  consideration  by  the 
Chief  of  the  Division  of  Materials  and  Tests  of  ISHC. 
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2 . 6  Desiccator 

A  desiccator,  about  ft  i richer,   (approximately  200  mm)  in  dia- 
meter containing  anhydrous  silica  pel  or  other  r.ui  table 
desiccant . 
3.  Sampling 

3.1  A  sample  weighing  about  SO  gm.  shall  be  taken  from  the  thor- 
oughly mixed  portion  of  the  material  passing  the  (No.  10)  2.00  mm 
sieve,  which  has  been  obtained  in  accordance  with  the  standard 
Method  of  Dry  Preparation  of  Disturbed  Hoil  Samples  for  Tests 
(AARHTO  T  87-721). 

3.2  The  sample  shall  be  free  from  contaminating  substances,  and 
thoroughly  representative  of  the  material  to  be  tested. 

h.    Procedure 

h.l     The  air-dried  soil  shall  be  placed  in  a  container  and  oven- 
dried  for  12  hours  or  to  constant  weight  at  110°C.  +  5°C. 
(230°F.  ±  9°F. ). 

k.2     A  representative  oven-dried  portion  of  5.00  gm. ±  0.005  gm. 
shall  be  placed  in  the  crucible  (two  crucibles  per  sample). 
The  crucibles  shall  then  be  placed  in  an  electric  muffle 
furnace  for  h   hours  at  1200°F.  ±  50°F. 

k.3     After  ignition  is  completed  at  1200°F.  the  samples  shall  be 
placed  in  a  desiccator  and  weighed  as  soon  as  cool. 
5.  Calculation  and  Report 

5.1  Report  the  organic  content  as  the  weight  per  cent  of  the  original 


sample  which  is  calculated  as  follows: 

T  n     T      ...  (B-A)    x  100 

Loss-On-Igmtion     = 

a 
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Where : 

A  =  Weight  of  sample  after  ignition.  (^.3) 
B  =  Weight  of  sample  originally  taken.  (U.2) 

A -2.   Method  of  Test  for  Determining  The  Slake  Durability  Index  of  Shale* 
Test  Method  No.  Ind.  503-77 

1.  Scope 

1.1  The  slake  durability  test  is  a  measure  of  the  abrasion  dura- 
bility of  a  shale  while  submerged  in  water.   This  weight  loss 
is  then  expressed  as  the  durability  number  ( I  ) . 

2.  Apparatus 

2.1  Oven 

A  thermostatically  controlled  drying  oven  capable  of  maintaining 
temperatures  of  110°C.  ±  5°C.   (230°F.  +  9°F.)  for  drying  shale 
samples . 

2.2  Balance 

A  balance  conforming  to  the  requirements  of  AASHTO  M231 ,  Class  D. 

2. 3  Containers 

Containers  for  moisture  content,  made  of  metal  or  other  suitable 
material,  with  close-fitting  lids  to  nrevent  loss  of  moisture 
prior  to  or  during  weighing. 

2 . h     Soak  Tank 

A  tank  of  sufficient  size  to  allow  complete  immersion,  of  the 
Apparatus  from  the  Slake  Durability  Test,  in  water. 

2 . 5  Slake  Durability  Apparatus 

The  Slake  Durability  Apparatus  shall  consist  of  a  motor  drive 
unit  mounted  on  a  base  board  and  capable  of  revolving  a  test 


*This  test  method  is  only  tentative  and  is  under  consideration  by  the  Chief 
of  the  Division  of  Materials  and  Tests  of  ISHC. 
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drum  at  a  speed  of  20  revolutions  per  minute.   The  test  drum 
shall  he  manufactured  of  brass  comprising  2.0  mm  standard  mesh 
cylinder  of  unobstructed  length  100  mm  and  iho   mm  diameter. 
Both  ends  of  the  cylinder  shall  be  solid  and  incorporate  suit- 
able driving  clogs.   One  side  plate  shall  incorporate  a  quick 
release  mechanism  to  permit  easy  sample  handling.   The  test 
drum  shall  be  supported  on  water  lubricated  bearings  allowing 
^0  mm  unobstructed  clearance  below  the  drum  and  a  trough  water 
level  20  mm  below  the  axis  of  the  drum. 
The  water  trough  shall  be  manufactured  by  Porspex 
3.  Test  Cample 

3.1  The  shale  material  shall  be  in  an  unweathered  condition. 

3.2  The  material  shall  he  at  nature  moisture  content. 

3.3  The  sample  shall  be  free  from  contaminating  substances,  and 
thoroughly  representative  of  the  material  to  be  tested. 

3.^   Forty  pieces  of  shale  shall  be  selected  such  that  they  are 

each  roughly  equi dimensional ,  and  weighing  between  60  gm.  and 

70  gm. 

k.    Procedure:   Slake  Durability  Prv,  (T,)d. 

d 

h.l     Two  samples  of  ten  representative  shale  pieces,  each  weighing 
600  gm.  to  TOO  gm.  ,  shall  be  placed  in  a  container  and  oven- 
dried  for  12  hours  or  to  constant  weight  at  110°C.  ±  5°r- 
(230°F.  +  9°F.)   (See  note  l).   The  moisture  shall  be  calculated 

and  recorded  as  (f!I,)d. 

d 

h.2     'Tie  shale  material  shall  be  allowed  to  cool  for  30  minutes  at 
room  temperature.   The  oven-dried  shale  material  then  shall  be 
placed  in  the  test  drum  and  the  weight  recorded  (see  note  l)  . 
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k.3     The  sample  and  test  drum  shall  be  immersed  in  the  water  hath 

of  the  Slake  Durability  Apparatus.   The  trough  water  level  shall 
be  two  {?)    cm.  below  the  axis  of  the  drum  (see  note  l). 

h.h     One  sample  shall  be  rotated  for  200  revolutions  and  the  second 
sample  for  500  revolutions . 

U.5  At  the  conclusion  of  section  h.h,   the  test  drum  and  sample  shall 
be  removed  from  the  water  bath  of  the  Slake  Durability  Apparatus 
and  oven-dried  for  12  hours  or  to  constant  weight  at  110  C. 
±   5°C  (230°F.  ±  9°F.),  and  the  weight  recorded  (see  note  l). 
5.  Procedure:   Soaked  Slaked  Durability  Index,  (I  )s. 

5.1  Two  samples  of  ten  representative  shale  pieces,  each  weighing 
600  gm.  to  700  gm. ,  shall  be  oven-dried  for  12  hours  or  to  con- 
stant weight  at  110°C.  ±  5°C   (230°F.  ±Q°F.  )   (see  note  l). 

The  moisture  shall  be  calculated  and  recorded  as  (MI,)s. 

d 

5.2  The  shale  material  shall  be  allowed  to  cool  for  30  minutes 

at  room  temperature.   The  oven-dried  shale  material  then  shall 
be  placed  in  the  test  drum  and  the  weight  recorded  (see  note  l). 

5.3  The  sample  and  test  drum  shall  be  immersed  in  the  Soak  Tank 
for  2h   hours.  At  the  conclusion  of  soaking,  the  sample  and 
test  drum  shall  be  immersed  in  the  water  bath  of  the  Slake 
Durability  Apparatus.   The  trough  water  level  shall  be  two  (2) 
cm.  below  the  axis  of  the  drum  (see  note  l). 

5.1+  One  sample  shall  be  rotated  for  200  revolutions  and  the  second 

sample  for  500  revolutions. 
5.5  At  the  conclusion  of  section  5.^,  the  test  drum  and  sample 

shall  be  removed  from  the  water  bath  of  the  Slake  Durabilitv 
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Apparatus  and  oven-dried  for  12  hours  or  to  constant  weight 
at  110°C.  +  5°C.  (230°F.  ±  9°F),  and  the  weight  recorded 
( see  note  1 )  . 
6.  Calculations  and  Report 

6.1  Report  the  moisture,  (MI  )d  or  (MI  )s,  as  the  weight  loss 
in  per  cent  of  the  oven-dried  material,  calculated  as 
follows : 

("ljd  or  (MI.)s,  Per  Cent  =  ^   -  „  x  100 

a  d  n    -   L 

Where:   A  =  Weight  of  container  and  natural  material 

B  =  Weight  of  container  and  oven-dried  material. 
C  =  Weight  of  container. 

6.2  Report  the  Slake  Durability  Index,  ( I  )d  or  (I  )s,  calculated 
as  follows: 

(ljd  or  (I  )s  =  ~~ ~-     x  100 
d        d       D  -  F 

Where:   D  =  Weight  of  drum  plus  oven-dried  material  before 

Test . 

K  =  Weight  of  drum  plus  oven-dried  material  retained 

after  test. 

F  =  Weight  of  clean  and  dry  drum. 

6.3  Report  the  classification  of  the  shale  material  with  respect 
to  Figure  A-l.  Durability  is  determined  based  on  the  500 
revolution  values.  The  200  revolution  values  shall  be  report- 
ed  as  additional  information. 


Note  1:   Care  should  be  taken  to  minimize  additional  breakdown  of  the 
material. 
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A- 3 .      Method  of  Test    for  Determining    the   Slaking    fndex   of  Shale 
Test   Method   No.    Ind.    502-73 

1.  Scope 

1.1  This  method  covers  the  determination  of  the  Slakinp;  Index 

to  he  used  in  classifying  shale.   The  procedure  is  based  on 
the  slaking  characteristics  of  shale  in  water  with  five  cycles 
of  wettinp;  and  drying. 

2.  Apparatus 

2.1  Oven 

A  thermostatically  controlled  drying  oven  capable  of  maintain- 
ing temperatures  of  110°C.  ±  5°C.  (230°F.  ±  9°F. )  for  dryinp 
shale  samples. 

2 . 2  Balance 

A  balance  conforming  to  the  requirements  of  AASHTO  M231,  Class  C. 

2 . 3  Beakers 

Pyrex  beakers  600  ml. 

3.  Test  Sample 

3.1  The  test  sample  shall  consist  of  six  pieces  of  shale  carefully 
selected  to  be  representative  of  the  shale  sample  submitted  for 
testing,  (Note  1. ) .   Each  piece  shall  be  approximately  equi -dimen- 
sional, shall  weigh  between  100  and  150  grams,  and  shall  be  free 
of  foreign  or  contaminating  substances.   Trimming  or  breaking  of 
pieces  to  obtain  approximately  equal  dimensions  is  permissible. 


Note  1.   Shale  samples  received  for  testing  shall  be  kept  in 
sealed  containers  to  retain  the  natural  moisture  con- 
tent to  minimize  weathering. 
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Procedure 

U.l   Each  piece  of  the  test  sample  shall  be  placed  in  a  separate 

beaker,  weighed,  then  oven-dried  to  constant  weight  at  110°C. 
+  5°C.   (230°F.  +  9°F.).   (Note  2).   The  moisture  content  shall 
be  calculated  and  recorded. 

h.2     The  shale  material  shall  be  allowed  to  cool  for  30  minutes  at 
room  temperature.   Distilled  water  then  shall  be  poured  into 
the  beakers  so  that  the  material  is  covered  by  at  least  1  cm. 
of  water.   (Note  3) . 

h.3     The  condition  of  the  shale  samples  shall  be  observed  at  5  minutes, 
30  minutes,  2  hours,  k   hours,  and  2k   hours  after  immersion.   The 
condition  of  the  material  shall  be  recorded  as :   complete  break- 
down, partial  breakdown,  or  no  change.   If  the  material  appears 
to  remain  intact,  the  cloudiness  of  the  water  shall  be  noted. 
If  there  is  variability  in  the  condition  of  the  6  test  pieces, 
the  variability  shall  be  recorded. 


Note  2.   Checking  every  moisture  content  sample  to  determine 

that  it  is  dried  to  a  constant  weight  is  impractical. 
In  most  cases,  drying  of  a  moisture  content  sample 
overnight  (15  or  16  hr.)  is  sufficient.   In  cases  where 
there  is  doubt  concerning  the  adequacy  of  overnight  dry- 
ing, drying  should  be  continued  until  the  weights  after 
two  successive  periods  of  drying  indicate  no  change  in 
weight . 

Note  3.   Care  shall  be  taken  to  minimize  additional  breakdown 
of  the  material. 
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h .k   At  the  conclusion  of  the  2h   hours  of  immersion  in  water,  the 
material  in  each  heaker  shall  he  washed  with  distilled  water 
on  a  2.00  mm  (No.  10)  sieve.   The  material  retained  on  the 
sieve  shall  then  he  washed  back  into  the  beaker,  decanted  and 
oven-dried  to  constant  weight  at  a  temperature  of  110°C.  ±  5°C. 
(230°F.  ±  9°F.)  and  weighed.   (Notes  2  and  3).  The  percent  loss 
from  initial  oven-dried  weight  shall  then  be  calculated  and 
recorded  (Slaking  index  for  that  cycle). 

k.3   Repeat  steps  U.2  through  h.k   for  five  complete  cycles  or  until 
the  material  is  completely  slaked,  whichever  occurs  first. 
Calculation  and  Report 

5.1  Calculate  the  moisture  content  as  follows: 

A-B 
Moisture,  percent  =       x  100 

B— C 

A  =  Total  weight  of  the  6  beakers  and  natural  material. 

B  =  Total  weight  of  the  6  beakers  and  oven-dried  material. 

C  =  Total  weight  of  the  6  beakers . 

Report  the  moisture  content  to  the  nearest  0.1%. 
5-2  Report  the  condition  of  shale  material  for  each  cycle  as  set 

out  in  k.3 
5.3  Calculate  the  Slaking  Index  at  the  end  of  each  cycle  as  follows: 

R— T) 
Slaking  Index  =  ■   -—  x  100 
B— L 

B  =  Total  Weight  of  the  6  beakers  and  oven-dried  material  of  k.l. 
C~=   Total  Weight  of  the  6  beakers. 

D  =  Total  Weight  of  the  6  beakers  and  oven-dried  material  retained 
on  the  2.00  mm  (No.  10)  sieve. 
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Report  the  Slaking  Index  for  each  cycle  to  the  nearest  0.1/T. 
5 .U   Report  the  Classification  of  the  shale  material,  with  respect 
to  Figure  A-l.   (Note  k) . 

A-U.   Method  of  test  for  determining  Fissility  Number  (Deo,  1972) 
1.   Breaking  Characteristics  of  Shales. 

The  breaking  characteristics  may  be  the  most  descriptive  feature  for 
shales.   They  can  be  classified  as  massive,  flaky-fissile  and  flaggy- 
fissile.   Fissility  is  associated  with  a  parallel  arrangement  of  clay 
particles,  and  non-fissility  with  a  random  arrangement.   The  nature  of 
cementing  agents  is  also  an  important  factor  in  influencing  fissility. 

Massive  rocks  have  no  preferred  directions  of  cleaving  and  breaking. 
Most  of  the  fragments  are  blocky.   Flaggy  rocks  will  split  into  fragments 
of  varying  thickness,  but  the  width  and  length  are  many  times  greater 
than  the  thickness,  and  the  two  essentially  flat  sides  are  approximately 
parallel.   Flaky  shales  split  along  irregular  surfaces  parallel  to  the 
bedding,  and  into  uneven  flakes,  thin  chips,  and  wedge-like  fragments 
whose  length  seldom  exceeds  three  inches. 

Shales  were  broken  by:   (a)  a  hammer  having  a  large  area  of  contact, 
and  (b)  by  striking  pieces  of  shale  against  each  other.   About  1000  gm 
of  shale  was  broken  in  this  way,  and  approximately  the  same  breaking 
effort  was  applied  to  each  shale. 

Shale  pieces  with  massive,  flaggy  and  flaky  characteristics  were 
visually  separated  and  weighed.   Proportions  of  the  three  different 


Note  h.     Additional  tests  are  necessary  to  classify  the  shale  material 
if  it  does  not  slake  completely. 
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breaking  types  were  determined  to  the  nearest  10  percent. 
2.  Determining  Fissility  Number. 

Flaky  and  flaggy  are  two  characteristics  of  fissility,  and  therefore 
the  fissility  number  should  be  some  weighted  sum  of  the  two.   Tt  was 
assumed  that  the  fissility  number  should  be  equal  to  the  percent  flaky 
component  plus  a  constant  times  the  percent  flaggy  component.   Though  the 
size  and  weight  of  flaggy  or  flaky  pieces  for  a  given  shale  varies  with 
the  breaking  effort,  the  flaggy  pieces  will  be  heavier  than  the  flaky 
pieces,  for  a  given  breaking  effort.   Typically  the  weight  of  flaky 
pieces  varied  between  5  and  100  percent  of  flaggy  pieces,  and  the  aver- 
age weight  of  flaky  pieces  was  0.35  times  the  average  weight  of  flaggy 
pieces . 

Therefore,  the  fissility  number  was  defined  as  the  sum  of  percent 
flakiness  and  0.35  times  percent  flagginess . 

A-5.  Method  of  test  for  determining  Modified  Soundness  (Deo,  1972). 

This  test  measures  the  degradation  of  shales  when  subjected  to  five 
cycles  of  alternate  wetting  and  drying  in  a  sodium  sulfate  solution. 

The  test  was  modified  from  AfiTM  C  88-63,  which  is  used  to  determine 
the  resistance  of  aggregates  to  disintegration  by  sodium  sulfate  or  mag- 
nesium sulfate.   The  standard  test  uses  a  fully  saturated  solution,  but 
this  is  too  severe  for  shales,  and  after  a  series  of  trials,  the  satur- 
ation was  reduced  to  50%. 

The  charge  of  shale  fragments  was  1000  gm,  of  which  330  gm  was  be- 
tween 1/2  in-   and  3/8  in.,   and  670  gm  was  between  3A  in.  and  1/2  in. 
Pieces  in  this  size  range  were  roughly  equidimensional.   Larger  pieces 
tended  to  be  plate  shaped,  due  to  the  laminated  nature  of  the  sediment. 


A-12 


Definition  of  size  by  a  sieving  process  of  course  becomes  more  arbitrary 
as  the  pieces  depart  from  a  bulky  shape.   The  sample  is  wa3h4d  with  water, 
and  oven  dried  at  105  to  110°C  before  weighing. 

A  saturated  solution  of  anhydrous  granular  sodium  sulfate  is  pre- 
pared in  accordance  with  AGTM  C  08-63  procedures.  The  solution  is  di- 
luted to  30%  saturation  by  adding  an  equal  ojnount  of  water.  The  solu- 
tion is  prepared  at  least  2h   hours  in  advance  of  the  start  of  test. 

The  sample  is  immersed  in  the  sodium  sulfate  solution  for  not  less 
than  16  hours  and  not  more  than  18  hours.   The  solution  covers   the  shale 
chunks  to  a  depth  of  at  least  1/2  in.   The  immersion  is  conducted  at  a  room 
temperature  of  72°  ±  2°F.   The  sample  is  removed  from  the  solution,  drained 
for  16  minutes,  placed  in  the  drying  oven  at  105  to  110°0 ,  and  dried  to 
constant  weight.   After  the  sample  has  cooled  to  room  temperature,  the 
process  is  repeated. 

Upon  completion  of  five  cycles  of  immersion  and  drying,  the  sample 
is  washed  with  water  until  free  of  sodium  sulfate,  as  determined  by  the 
reaction  of  the  wash  water  with  barium  chloride  (BaCl  ).   It  is  then  dried 
and  fractioned  on  a  5/l6  in.  sieve.   The  weight  retained  on  the  sieve  is 
determined.   Each  test  is  repeated  at  least  once,  and  average  values  are 
reported. 

The  Soundness  Index,  is  defined  as  the  percent  retained  by  weight  on 

the  5/l6  in.  sieve.   Durability  is  considered  to  increase  with  increase 

in  I  value, 
s 
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DATA  SEGREGATED  ACCORDING  TO  GEOLOGICAL  FORMATION? 
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